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The electromagnetic properties of the baryon octet are calculated in quenched QCD on a 20 3 x 
40 lattice with a lattice spacing of 0.128 fm using the fat-link irrelevant clover (FLIC) fermion 
action. FLIC fermions enable simulations to be performed efficiently at quark masses as low as 
300 MeV. By combining FLIC fermions with an improved-conserved vector current, we ensure 
that discretisation errors occur only at 0(a 2 ) while maintaining current conservation. Magnetic 
moments and electric and magnetic radii are extracted from the electric and magnetic form factors 
for each individual quark sector. From these, the corresponding baryon properties are constructed. 
Our results are compared extensively with the predictions of quenched chiral perturbation theory. 
We detect substantial curvature and environment sensitivity of the quark contributions to electric 
charge radii and magnetic moments in the low quark mass region. Furthermore, our quenched 
QCD simulation results are in accord with the leading non-analytic behaviour of quenched chiral 
perturbation theory, suggesting that the sum of higher-order terms makes only a small contribution 
to chiral curvature. 
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I. INTRODUCTION 

The study of the electromagnetic properties of baryons 
provides valuable insight into the non-perturbative struc- 
ture of QCD. Baryon charge radii and magnetic moments 
provide an excellent opportunity to observe the chiral 
non-analytic behaviour of QCD. Although the first cal- 
culations of hadronic electromagnetic form factors ap- 
peared almost 20 years ago 0, y, 13 > until recently the 
state-of-the-art calculations of the electromagnetic prop- 
erties of octet 0, |5j and decuplet [6j baryons and their 
electromagnetic transitions [?J appeared almost 15 years 
ago. 

However, over the last couple years there has been 
an increase in activity in the area of octet baryon elec- 
tromagnetic structure, mainly by the Adelaide group 
H ©, Il3jp and the QCDSF [l| and LHPC collabo- 
rations (l3j |. In this paper we improve upon our prelim- 
inary results reported in Ref. [8j and describe in detail 
the origin of the lattice simulation results featured in 
Refs. 0, E3 an d UM determining the strangeness mag- 
netic moment and charge radius of the nucleon respec- 
tively. 

The extraction of baryon masses and electromagnetic 
form factors proceeds through the calculation of Eu- 
clidean two and three-point correlation functions, which 
are discussed at the hadronic level in Section Til Bl and 
at the quark level in Sections III CI and III Dl Throughout 
this analysis we employ the lattice techniques introduced 
in |j| . We briefly outline the main aspects of these tech- 
niques in section ITTT1 The correlation functions directly 
proportional to the electromagnetic form factors of inter- 
est are analysed in Sec. II VI The results are presented and 



discussed in Section where an extensive comparison is 
made with the predictions of quenched chiral perturba- 
tion theory (QxPT) [lj, [l5| . A summary and discussion 
of future work is provided in Section IVll 



II. THEORETICAL FORMALISM 

A. Interpolating fields 

In this analysis we work with the standard established 
interpolating fields commonly used in lattice QCD simu- 
lations. The notation adopted is similar to that of |4|. To 
access the proton we use the positive parity interpolating 
field 



X p+ (x) 



r abc („,aT 



(u aT (x) C l5 d b (x)) u c (x) , (2.1) 



where the fields u, d are evaluated at Euclidean space- 
time point x, C is the charge conjugation matrix, a, b 
and c are colour labels, and the superscript T denotes 
the transpose. In this paper we follow the notation of 
Sakurai. The Dirac 7 matrices are Hermitian and satisfy 

{7/*>7«/} = 2( W witu en* = ( 1 / 2 *)[7p)7^]- This interpo- 
lating field transforms as a spinor under a parity trans- 
formation. That is, if the quark fields q a {x) (q = u, d, ■ ■ •) 
transform as 

Vq a (x)V^ =+ 7 o9 a (£), (2.2) 

where x = (xq, —x), then 

V X p+ (x)V^+ 7 oX p+ (x). (2.3) 

The neutron interpolating field is obtained via the ex- 
change u «-► d, and the strangeness —2, 5 interpolating 



fields are obtained by replacing the doubly represented u 
or d quark fields in Eq. 12. If) by s. Similarly, the charged 
strangeness — 1, S interpolating fields are obtained by 
replacing the singly represented u or d quark fields in 
Eq. 12.1f> by s. For the S° hyperon one uses |4| 

X S °(x) = i=e o6c { (u aT (x) C l5 s b (x)) d c (x) 



(d aT {x)C lb s b {x))u c {x))(2A) 



Note that x s transforms as a triplet under SU(2) 
isospin. An SU(2) isosinglet interpolating field for the 
A can be constructed by replacing " + " — > " — " in 
Eq. (|21}. For the SU(3) octet A interpolating field, one 
has 



XV) = -^=e ahc {2 (u aT {x) C l5 d b {x)) s c {x) + (u aT {x) C l5 s b (x)) d c {x) - (d aT (x) C l5 s b {x)) u c (x)} . (2.5) 



We select this interpolating field for studying the A in the following. 



B. Correlation functions at the hadronic level 



The extraction of baryon masses and electromagnetic form factors proceeds through the calculation of the ensemble 
average (denoted (••■)) of two and three-point correlation functions. The two-point function is defined as 



< G BB (t;p,T) )=J2 e ~*' 3rf>a < Q \T( X a (x)x (O)) | Q ) 



(2.6) 



Here $7 represents the QCD vacuum, T is a 4 x 4 matrix in Dirac space and a, (3 are Dirac indices. At the hadronic 
level we insert a complete set of states B,p,s ) and define 



< fi | X (0) | B,p, S )=Z B (p)J^-u(p,s) , 



(2.7) 



where Zb(p) represents the coupling strength of x(0) to baryon B, and E p = \JpP- + M 2 . A momentum dependence 
for Zb(p) is included for the case where a smeared sink is employed. For large Euclidean time 



(G BB (t;p,T)) 



Z B (p)Z B (p) ^ Et ,t 



2E„ 



e -& P t tr [r(-fy • p + M)] 



(2.8) 



Here Zb(p) is the coupling strength of the source x(0) to the baryon. Again, the momentum dependence allows for 
the use of smeared fermion sources in the creation of the quark propagators and the differentiation between source 
and sink allows for our use of smeared sources and point sinks in the following. Similarly the three-point correlation 
function for the electromagnetic current, j'^(x), is defined as 



< G B ^ B (t 2 ,t i; p',p;r) >= Y, e-^e+'^Hr^ < fi | T (x a (x 2 )j' i (x 1 )x (O)) \ fJ ) 



(2.9) 



For large Euclidean time separations ti — t\ » 1 and t\ » 1, the three-point function at the hadronic level is 
dominated by the contribution from the ground state 

(G B ^ B {t 2 ,t 1 - ) p',p;T))=Y e ~ Ep ' {t2 ~ tl)e ~ EptlTl3a {n\x a \p',s')(p',s' \j" \p,s)(p,s\x 13 |fi) . (2.10) 
The matrix element of the electromagnetic current has the general form 

/ ' ' I -M I \ 

\p ,s \r \p,s) = 



M 2 N ' 2 

E p E p / 



,2\„ui/ 1 



u(p',s') ( FrteV " F ^)^m I U(P ' S) 



(2.H) 



where q = p' — p. To eliminate the time dependence of the three-point functions we construct the following ratio, 

1/2 



R(t2,ti;p',p;T,r';^) = 



( G B ^ B (t 2 ,t i;P ',p;T) )( G B ^ B (t 2 ,t i; ~p,-pi;T) ) 



(2.12) 



(G BB (t 2 ;p'-X))(G BB (t 2 :-p;T')) 

For large time separations t 2 — t\ » 1 and t\ » 1 this ratio is constant in time and is proportional to the 
electromagnetic form factors of interest. We further define a reduced ratio R(j>' ,p]T,T ; fj.) as 

1/2 



#(p\p;i\r -p.) = 



2E n 



E p + M 



1/2 



2E,„ 



E.. 



M 



R(t2,ti;p',p;T,T ;/i) , 



(2.13) 



from which the Sachs forms for the electromagnetic form 
factors 



r 4 = - 



i / / o 



o o 



(2.19) 



G E (q 2 



Fi(Q 2 ) 



(2M) 2 
G M (q 2 ) = F^q 2 ) + F 2 (q 2 ) , 



F 2 (q 2 



(2.14) 
(2.15) 



may be extracted through an appropriate choice of T and 
r . A straight forward calculation reveals 

g E (q 2 ) = l(g,o;r 4 ,r 4 ,4), (2.16) 

\e l3k q l \GM(q 2 ) = {E q + M)R(q, 0;^-, T 4 , k) , (2.17) 

\q k \Q E {q 2 ) = (^ + Af)i?(g,0;r 4 ,r 4 ,fc),(2.18) 



where 



3 ~ 2 I 



C. Correlation functions at the quark level 



Here the two and three-point functions of Sec. Ill Bl are 
calculated at the quark level by using the explicit forms 
of the interpolating fields of Sec. Ill Al and contracting out 
all possible pairs of quark field operators. These become 
quark propagators in the ensemble average. For conve- 
nience, we introduce the shorthand notation for the cor- 
relation functions Q of quark propagators S 



G{ s fi> s f2> s h) = e<1 Ce ° 



S a f f (x, 0) tr S% T (x, 0)S%' (x, 0) + S a f f (x, 0) S% T (x, 0) S%' (x, 0) 



(2.20) 



where S 1 ?" (x, 0) are the quark propagators in the background link-field configuration U corresponding to flavours 
/i_3. This allows us to express the correlation functions in a compact form. The associated correlation function for 
X p+ can be written as 



G p+ {t,p;T) = /^V l P*tr 



r G ( s u , cSdC , s u 



(2.21) 



where (• • •) is the ensemble average over the link fields, T is the T± projection operator that separates the positive 
and negative parity states, and C = C75. For ease of notation, we will drop the angled brackets, (•••), and all the 
following correlation functions will be understood to be ensemble averages. 

Two-point correlation functions for other octet baryons composed of a doubly-represented quark flavour and a 
singly-represented quark flavour follow from Eq. (|2.21l) with the appropriate substitution of flavour subscripts. The 
correlation function for the neutral member S° is given by the average of correlation functions for the charged states 
E + and £~. Finally the correlation function for A obtained from the octet-interpolating field of Eq. (|2.5ll is 

G AS (t,p; r) = i £ e-<**tr r j 2 Q (s„ CS u C-\ S d ) + 2 Q (s„ CS d C~\ S u 



+ 2g(s d ,cs u c- 1 ,s s \ 2 

— G [ Sd, CS S C~ , S u 



GbjG 1 Of 



G (S u , 

g{s a) cs a c- x , s d y> 



(2.22) 



\u,d,s 



FIG. 1: Diagrams illustrating the two topologically different 
insertions of the current within the framework of lattice QCD. 



D. Three-point functions at the quark level 

In determining the three point function, one encoun- 
ters two topologically different ways of performing the 
current insertion. Figure ^ displays skeleton diagrams 
for these two insertions. These diagrams may be dressed 



J 



with an arbitrary number of gluons (and additional sea- 
quark loops in full QCD). Diagram (a) illustrates the 
connected insertion of the current to one of the quarks 
created via the baryon interpolating field. This simple 
skeleton diagram does indeed contain a sea-quark compo- 
nent, as upon dressing the diagram with gluon exchange, 
quark-loop and Z-diagrams flows become possible. It is 
here that "Pauli-blocking" in the sea contributions, cen- 
tral to obtaining violation of the Gottfried sum rule, are 
taken into account. Diagram (b) accounts for an alterna- 
tive quark-field contraction where the current first pro- 
duces a disconnected q q loop-pair which in turn interacts 
with the valence quarks of the baryon via gluons. 

Thus, the number of terms in the three-point function 
is four times that in Eq. (|2.21(l . The correlation function 
for proton matrix elements obtained from the interpola- 
tor of Eq. H2.1|) is 



q(s u (x 2 ,x 1 ,0), CS d (x 2 ,0)C-\ S u (x2,0))+g(s u (x 2 ,0), CS d (x 2 ,0)C-\ S u (x 2 ,x u 0) 
+g(s u {x 2 ,0), CS d {x 2 ,x 1 ,0)C'\ S u (x 2 ,0) 



J2 e q J2tr[S i j(x 1 ,xi) 7fi ]g[S a (x2,0),CS d (x 2 ,0)C- 1 , S u (x 2 ,0)j , 

q=u, d, s i 



(2.23) 



r 



where 

S a q a '(x 2 ,x 1 ,0)^e q J2s^(x 2 ,x 1 )^Sl a \x 1 ,0), 

(2.24) 
denotes the connected insertion of the electromagnetic 
current to a quark of charge e q . 

The first two terms of Eq. (|2.23|) provide the connected 
insertion contribution of the u-quark sector to the pro- 
ton's electromagnetic properties, whereas the third term 
provides the connected d-quark contribution. The latter 
term of Eq. (I2.23() accounts for the "disconnected" loop 
contribution depicted in Fig. \T]p. Here, the sum over the 
quarks running around the loop has been restricted to 
the flavours relevant to the ground state baryon octet. 
In the SU(3)-fiavour limit the sum vanishes for the elec- 
tromagnetic current. However, the heavier strange quark 
mass allows for a non-trivial result. 

The "disconnected" current insertion requires a numer- 
ical estimate of S™(xi,Xi) for the lattice volume of di- 
agonal spatial indices. As this requires numerous source 
vectors in the fermion-matrix inversion, determination 
of this propagator is numerically intensive [l|| [13, 13 • 
Indeed, an indirect method using experimental results, 
chiral effective field theory and the lattice results from 
the connected current insertion presented herein, pro- 
vides the most precise determinations of these quark loop 



contributions to the nucleon's electromagnetic structure 
UaUIlLUl & t present. This approach should be viewed as 
complementary to an ab initio determination via lattice 
QCD which awaits a next-generation dynamical-fermion 
simulation of QCD [l4| . 

It is interesting to examine the structure of the con- 
nected insertion contributions to the proton's structure. 
Here, we see very different roles played by u and d quarks 
in the correlation function, in that only the d-quark ap- 
pears in the second position of Q. The absence of equiv- 
alence for u and d contributions allows the connected 
quark sector to give rise to a non-trivial neutron charge 
radius, a large neutron magnetic moment, or a violation 
of the Gottfried sum rule. As each term of Eq. I|2.23|l can 
be calculated individually, it is a simple task to isolate 
the quark sector contributions to the baryon electromag- 
netic properties. 

Another interesting point to emphasise, is that there is 
no simple relationship between the properties of a partic- 
ular quark flavour bound in different baryons. For exam- 
ple, the correlator for E + is given by l|2.23|l with d—>s. 
Hence, a u-quark propagator in E + is multiplied by an 
s-quark propagator, whereas in the proton the u-quark 
propagators are multiplied by a d-quark propagator. The 
different mass of the neighbouring quark gives rise to an 
environment sensitivity in the u-quark contributions to 



obser vables 0, IE 0, S E3, Eol EH ■ This point sharply con- 
trasts the concept of an intrinsic quark property which is 
independent of the quark's environment. This concept of 
an intrinsic quark property is a fundamental foundation 
of many constituent based quark models and is not in 
accord with QCD. 



III. LATTICE TECHNIQUES 

A. Gauge and quark actions 

The simulations are performed using the mean-field 
O(a 2 )-improved Luscher-Weisz 22] plaquette plus rect- 
angle gauge action on a 20 3 x 40 lattice with periodic 
boundary conditions. The lattice spacing a = 0.128 fm 
is determined by the Sommer scale ro = 0.50 fm |23j . 
This large volume lattice ensures a good density of low- 
lying momenta which are key to giving rise to chiral non- 
analytic behaviour in the observables simulated on the 

lattice HE 03. 

We perform a high-statistics analysis using a large 
sample of 400 configurations for our lightest eight quark 
masses. We also consider a subset of 200 configurations 
for our three heaviest quark masses to explore the ap- 
proach to the heavy-quark regime. A sub-ensemble bias 
correction is applied multiplicatively to the heavy quark 
results, by matching the central values of the 200 con- 
figuration sub-ensemble and 400 configuration ensemble 
averages at n = 0.12780. The error analysis is performed 
by a third-order, single-elimination jackknifc. 

For the quark fields, we use the Fat-Link Irrelevant 
Clover fermion action 24] 



cFL 



cFL 



i g C S w k r 

2« L ) 4 



ip{x)a IJ , v F tlv 'ip(x), (3.1) 



where F^ v is an O(a 4 )-improved lattice definition J2|| 
constructed using fat links and itg L is the plaquette mea- 
sure of the mean link calculated with fat links. The mean- 
field improved Fat-Link Irrelevant Wilson action is 



Sw = Ex ^(XMX)+K^(X) 



7m — WC + aO 

V "0 



UUx - fi) 



wo 
U^(x-fi) 



ip(x-fi))-r[ % L ip(x + (i) 



,FL 



ip(x- A) 



(3.2) 



with k = l/(2m + 8r). We take the standard value r = 1. 
Our notation uses the Pauli representation of the Dirac 
7-matrices |26| . where the 7-matrices are hermitian and 
C/iw = [7^ 1 7i/]/(2i). Fat links are constructed by per- 
forming tt-ape = 6 sweeps of APE smearing, where in 
each sweep the weights given to the original link and the 
six transverse staples are 0.3 and (0.7/6) respectively. 
The FLIC action is closely related to the mean-field im- 
proved clover (MFIC) fermion action in that the latter 



is described by Eqs. (|3.1(l and i|3.2|) with all fat-links re- 
placed by untouched thin links and F^ v defined by the 
1 x 1-loop clover definition. 

For fat links, the mean link u ~ 1, indicating that 
perturbative renormalisations are small for smeared links 
and are accurately accounted for by small mean-field im- 
provement corrections. As a result, mean-field improve- 
ment of the coefficients of the clover and Wilson terms of 
the fermion action is sufficient to accurately match these 
terms and eliminate O(a) errors from the fermion action 
|27j . An added advantage is that access to the light quark 
mass regime is enabled by the improved chiral properties 
of the FLIC fermion action [28J] • 

Time slices are labeled from 1 to 40, and a fixed bound- 
ary condition at t = 40 is used for the fermions. An anal- 
ysis of the pion correlator indicates that the effects of this 
boundary condition are negligible for t < 30, and all of 
our correlation-function fits are performed well within 
this regime. 

Gauge- invariant Gaussian smearing [23, y(| in the spa- 
tial dimensions is applied at the source at t = 8 to in- 
crease the overlap of the interpolating operators with the 
ground state while suppressing excited state contribu- 
tions. 

Table [I] provides the kappa values used in our simula- 
tions, together with the calculated 7r and octet baryon 
masses. While we refer to m 2 . to infer the quark masses, 
we note that the critical value where the pion mass van- 
ishes is K cr = 0.13135. 

We select k — 0.12885 to represent the strange quark 
in this simulation. At this k the ss pseudoscalar mass is 
0.697 GeV, which compares well with the experimental 



value of 2 j 



(0.693 GeV) 2 , motivated by leading 



order chiral perturbation theory. 



B. Improved conserved vector current 

For the construction of the C(a)-improved conserved 
vector current, we follow the technique proposed by Mar- 
tinelli et al. |3l| . The standard conserved vector current 
for Wilson-type fermions is derived via the Noether pro- 
cedure 



ft 



= j[#c)(7m - r)U,,(x)ilj(x + p.) 

+ ^(x + fi)( 1(1 +r)Ul{x)^(x) 
+ (x -* X - /})] . 



(3.3) 



The C(a)-improvemcnt term is also derived from the 
fermion action and is constructed in the form of a to- 
tal four-divergence, preserving charge conservation. The 
O(a)-improved conserved vector current is 



Jfi — 



3°{x) + r -CcvcaJ2dp^( x )°pMx)) , (3-4) 



TABLE I: Hadron masses in appropriate powers of GeV for various values of the hopping parameter, n. 
experimentally measured values are indicated at the end of the table. 



For reference, 



K 


ml 


N 


A 


£ 


B 


0.12630 


0.9972(55) 


1.829(8) 


1.728(10) 


1.700(9) 


1.612(11) 


0.12680 


0.8947(54) 


1.763(9) 


1.681(10) 


1.656(10) 


1.586(12) 


0.12730 


0.7931(53) 


1.695(9) 


1.632(11) 


1.566(11) 


1.558(12) 


0.12780 


0.6910(35) 


1.629(10) 


1.584(10) 


1.570(10) 


1.531(10) 


0.12830 


0.5925(33) 


1.554(10) 


1.530(10) 


1.521(10) 


1.502(10) 


0.12885 


0.4854(31) 


1.468(11) 


1.468(11) 


1.468(11) 


1.468(11) 


0.12940 


0.3795(31) 


1.383(11) 


1.406(11) 


1.417(11) 


1.435(11) 


0.12990 


0.2839(33) 


1.301(11) 


1.347(11) 


1.371(11) 


1.404(11) 


0.13025 


0.2153(35) 


1.243(12) 


1.303(12) 


1.341(12) 


1.382(11) 


0.13060 


0.1384(43) 


1.190(15) 


1.256(13) 


1.313(12) 


1.359(11) 


0.13080 


0.0939(44) 


1.159(21) 


1.226(16) 


1.296(14) 


1.346(11) 


experiment 


0.0196 


0.939 


1.116 


1.189 


1.315 



where Ccvc is the improvement coefficient for the con- 
served vector current and we define 

d p (${x)^{x)) = ${x){V p + ^p)iP(x) , (3.5) 

where the forward and backward derivatives are defined 

as 



iKs)^* = ^(x + p)Ul(x) 



Ufa- fi)ip(x- A)] 

tyufa) 

ip(x - fi) U^x - A)] 



The terms proportional to the Wilson parameter r in 
Eq. I|3.3|) and the four-divergence in Eq. (|3.4(1 have their 
origin in the irrelevant operators of the fermion action 
and vanish in the continuum limit. Non-perturbative im- 
provement is achieved by constructing these terms with 
fat-links. As we have stated, perturbative corrections are 
small for fat-links and the use of the tree-level value for 
Ccvc = 1 together with small mean-field improvement 
corrections ensures that 0(a) artifacts are accurately re- 
moved from the vector current. This is only possible 
when the current is constructed with fat-links. Other- 
wise, Ccvc needs to be appropriately tuned to ensure 
all 0(a) artifacts are removed. 

In order to suppress contributions from excited states, 
large Euclidean times are required, both following the 
source at to, and following the current insertion at t\. 
Our two-point function analysis indicates that the ground 
state is isolated well by t = 14, largely due to an excellent 
selection for the source smearing parameters. Therefore 
the current insertion is performed at t\ = 14. 



C. Improved unbiased estimators 

The two and three-point correlation functions are de- 
fined as averages over an infinite ensemble of equilibrium 
gauge field configurations, but are approximated by an 
average over a finite number of configurations. To min- 
imise the noise in the results, we exploit the parity of the 
correlation functions |32| 

G(P,p, q: T) = s P G(-p', -p, -q; T), a P = ±l, 

J3.6) 
and calculate them for both p, p', q and — p, —p', —q. 
While this requires an extra matrix inversion to deter- 
mine S(x2,0;ti, — q, fi), the ratio of three- to two-point 
functions is determined with a substantial reduction in 
the statistical uncertainties. The improvement is better 
than that obtained by doubling the number of configura- 
tions. 

Similarly, the link variables {U} and {U*} are gauge 
field configurations of equal weight, and therefore we ac- 
count for both sets of configurations in calculating the 
correlation functions |33j . With the fermion matrix prop- 
erty 



M({U*})= [CM({U})C 



(3.7) 



it follows that 



S(x, 0; {U*}) = (dS(x, 0; {[/DC" 1 ) * , (3.8) 



S(x, 0; t, q, ^, {U*}) = [CS(x, 0; t, -q, /i; {U})C- 

(3.9) 
and therefore the correlation functions are purely real 
provided 



T = sc[CTC 



and sc — sp. 



(3.10) 



These conditions are satisfied with the selections for T 
indicated in Eq. (|2.19l) . In summary, the inclusion of 
both {[/} and {U*} configurations in the calculation of 
the correlation functions provides an improved unbiased 
estimate of the ensemble average properties incorporat- 
ing parity symmetry and significantly reducing statistical 
fluctuations. 



D. Fit regime selection criteria 

In fitting the correlation functions, the covariance- 
matrix based x 2 P er degree of freedom (dof) plays a 
central role. 

The correlated x 2 /dof is given by 



X 

dof N t - M 



, N t N t 

—T.T. (md 

(y(U) ~ T(U)) C-\u,ti) {y{t 3 ) - T{t 3 )) , 



where, M is the number of parameters to be fitted, Nt is 
the number of time slices considered, y(U) is the config- 
uration average value of the dependent variable at time 
U that is being fitted to a theoretical value T(ti), and 
C(ti,tj) is the covariance matrix. The elements of the 
covariance matrix are estimated via the jackknife method 



G (t i , tj 



N r -1 






N c ^ 

771 — 1 



(3.12) 



[y^(U) - y(ti)] [y^(tj) - y{tj)] > 
= (N c -l)x (3.13) 

tj- Yl ym{u) y^tti) - yfc) yfc) \ 



N r 



where N c is the total number of configurations and y^{t{) 
is the jackknife ensemble average of the system after re- 
moving the mth configuration. y(£j) is the average of all 
such jackknife averages, given by 



y(ti 



1 N ^ 

jj- c Y, y™^ 



(3.14) 



In the process of selecting the fit regimes, numerous fits 
are performed over a variety of start times and a variety 
of time durations. The following criteria are taken into 
account in selecting the preferred fit regime: 

1. The x 2 /dof is monitored and plays a significant 
role in determining the start time of the fit. Values 
within the range 0.5 to 1.5 are preferred and it is 
often possible to select a regime providing a perfect 
fit measure of 1. Start times for which the x 2 /dof 
increases significantly as the duration of the regime 
is increased are discarded. In practice, the x 2 /dof 
sets a lower bound for the start time, and other 
criteria may lead to a later start time for the fit. 



2. In some cases a monotonic systematic drift can 
be observed in the ratio of correlation functions 
H2.12J) which otherwise should be constant in time. 
Often the drift is sufficiently small to provide a 
X 2 /dof < 1.5. In these few cases, a later start 
time is selected to ensure that sufficient Euclidean 
time evolution has occurred to accurately isolate 
the ground state, suppressing systematic errors at 
the expense of larger statistical errors. 

3. As the quarks become lighter, the spacing between 
the ground and first excited states of the baryon 
spectrum becomes larger [23, lLJ , due to the more 
rapid reduction in the mass of the lower-lying state. 
This provides improved exponential suppression of 
excited-state contaminations. Hence, as one ap- 
proaches the light quark mass regime, a monotonic 
reduction in the starting time-slice of the fit regime 
may be possible. 

4. As the quark masses become lighter, the signal is 
lost to noise at earlier times. Hence the final time 
slice of the fit window is also monotonically de- 
creased as the quarks become lighter. We typically 
consider fit regimes of three to five time slices and 
preferably the latter when the signal is not obvi- 
ously lost to noise. 

5. For quark masses lighter than the strange quark 
mass, the splittings between adjacent quark masses 
are calculated and fit using the same techniques. 
By considering adjacent splittings, excited-state 
contributions which are less dependent on the 
quark mass (see item 3 above) are suppressed. In 
practice, good x 2 /dof are found one to two time 
slices earlier. 



A precise examination of the environment sensitivity 
of quark-sector contributions to baryon electromagnetic 
properties lies at the core of this investigation. For ex- 
ample, the doubly-represented u quark in the proton is 
to be compared with the doubly-represented u quark in 
S + ; the singly- represented u quark in the neutron with 
the u quark in 2°. Similarly, it is interesting to compare 
the strange and light quark sectors of S~ with those of A. 
Conventional models reverse the ordering of the observed 
magnetic moments. After the consideration of the pre- 
ceding criteria, the fit regimes are unified for each of the 
quark sector contributions wherever possible. This com- 
parison is done for each value of n governing the quark 
mass, reducing systematic errors associated with choos- 
ing different time-fitting regimes for similar quantities. 
For example, for the case of the doubly-represented u 
quark in the proton and E + , it is possible to equate the 
fit regimes for all but the two lightest quark masses where 
the x 2 1 dof insists on different fit regimes. 
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FIG. 2: Masses of the octet baryons. The SU(3) flavour limit 
is evident. Points on the y-axis indicate the experimentally 
measured masses for reference. 



IV. CORRELATION FUNCTION ANALYSIS 

The following calculations are performed with p = 
0, p' = q = \q\ x at q x a = 2t:/L x with L x = 20, 
the minimum non-zero momentum available on our lat- 
tice. We introduce Q 2 — —q 2 , as q 2 is negative (space- 
like). While Q 2 is dependent on the mass of the baryon, 
we find this mass dependence to be small. Indeed all 
form factors may be regarded as being calculated at 
Q 2 = 0.227±0.002 GeV 2 where the error is dominated by 
the mass dependence of the target baryon. Where a spa- 
tial direction of the electromagnetic current is required, 
it is chosen to be the z-dircction. Electric and magnetic 
form factors are extracted from our correlation functions 
as described in Eqs. (|2.16|) and IJ2.17JI . 



A. Baryon masses 

The masses of the baryon octet are plotted against 
m 2 in Fig. |2 and are tabulated in Table U We observe 
the SU(3)f limit at our sixth quark mass. The mass 
splitting between £ and A at the lowest pion mass (m^ = 
0.3064± 0.0072 GeV) on our lattice is 69±2 MeV which is 
only slightly smaller than the experimentally measured 
splitting of 76 MeV. Hence the generic features of the 
baryon-octet mass spectrum is reproduced well in our 
quenched simulation. 



B. Form factor correlators 

In general, the baryon form factors are calculated on a 
quark-sector by quark-sector basis with each sector nor- 
malised to the contribution of a single quark with unit 
charge. Hence to calculate the corresponding baryon 
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FIG. 3: Electric form factor of the proton and its quark sec- 
tors (including charge and quark number factors) at Q 2 = 
0.227(2) GeV 2 as a function of Euclidean time (£2) for m 2 = 
0.4854 GeV 2 , the SC/(3)-fiavour limit. The lines indicate the 
fitting windows and the best fit value. 



property, each quark sector contribution should be multi- 
plied by the appropriate charge and quark number. Un- 
der such a scheme for a generic form factor /, the proton 
form factor, / p , is obtained from the u- and d-quark sec- 
tors normalised for a single quark of unit charge via 

/ P = 2x|x/„ + lx (-^\ xf d . (4.1) 

The electric form factor of the proton and contribu- 
tions from the u- and d-quark sectors are plotted in Fig.EI 
as a function of Euclidean time at the 5J7(3)-flavour 
limit. Here, charge and quark number factors have been 
included such that the proton result is simply the sum 
of the illustrated quark sectors. The lines indicate the 
time slices selected for the fit using the considerations of 
Sec. l77TDl 

We find that substantial Euclidean time evolution is 
required following the current insertion to obtain accept- 
able values of the x 2 /dof; in this case seven time slices 
following the current insertion at t\ = 14. 

For light quark masses lighter than the strange quark 
mass, we fit the change in the form factor ratios of 
Eq. I|2.12|) from one quark mass to the next and add 
this to the previous result at the heavier quark mass. 
Figure 0] shows the quark sector contributions (including 
charge and quark number factors) to the electric form 
factor of the proton at Q 2 = 0.227(2) (GeV 2 ) as a func- 
tion of Euclidean time, £2 , for the ninth quark mass where 
to 2 = 0.2153(35) GeV 2 . The correlator is obtained from 
the splitting between the ninth and eighth quark mass 
states. The improvement of the plateau is apparent in 
Fig. 0] Still substantial Euclidean time evolution is re- 
quired to obtain an acceptable \ 2 /dof '. The onset of 
noise at this lighter quark mass is particularly appar- 
ent at time slice 24 for the d sector. Tables ITT1 to fVl list 
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FIG. 4: Quark sector contributions (including charge and 
quark number factors) to the electric form factor of the proton 
at Q 2 = 0.227(2) (GeV 2 ) as a function of Euclidean time, £2, 
for the ninth quark mass where m 2 = 0.2153(35) GeV 2 . The 
correlator is obtained from the splitting between the ninth 
and eighth quark mass states. The lines indicate the fitting 
windows and the best fit value. 
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FIG. 5: Magnetic form factors of H and its quark sec- 
tors (including charge and quark number factors) at Q 2 = 
0.227(2) GeV 2 as a function of Euclidean time (£2) for m 2 = 
0.4854 GeV 2 , the S'f7(3)-flavour limit. The lines indicate the 
fitting windows and the best fit value. 
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FIG. 6: Quark sector contributions (including charge and 
quark number factors) to the magnetic form factor of H° at 
Q 2 — 0.227(2) (GeV 2 ) as a function of Euclidean time, £ 2 , 
for the ninth quark mass where m 2 = 0.2153(35) GeV 2 . The 
correlator is obtained from the splitting between the ninth 
and eighth quark mass states. The lines indicate the fitting 
windows and the best fit value. 



the physical nucleon mass appears, has been done by 
multiplying the lattice form factor results by the ratio 
misr/ms- In this way the form factors are presented in 
terms of a constant unit; i.e. the nuclear magneton. 

The negative contribution of the u quark to the total 
form factor indicates that its spin is on average opposite 
to that of the doubly represented s quarks. This, as well 
the relative magnitude of the contributions, is in quali- 
tative agreement with simple constituent quark models 
based on SU(6) spin-flavour symmetry. 

In Fig. we present the Euclidean time dependence 
of the the magnetic form factors of S° calculated at the 
ninth quark mass where to 2 . = 0.2153(35) GeV 2 . Again 
the early onset of acceptable plateau behaviour is appar- 
ent here. 

Results for the quark-sector contributions to the mag- 
netic form factors of octet baryons are summarised in 
Tables EH to IDO 



DISCUSSION OF RESULTS 



the electric form factors for all the octet baryons at the 
quark level for the eleven quark masses considered. In 
the tables, the selected time frame, the fit value and the 
associated X 2 /dof are indicated. 

Turning now to the magnetic form factors, Fig.[S]shows 
the magnetic form factor of S° and its quark sectors (in- 
cluding charge and quark number factors) as a function 
of Euclidean time at the 5C/(3)-fiavour limit. Preferred 
fit windows following from the criteria of Sec. IIII Dl and 
best fit values are indicated. 

Here the conversion from the natural magneton, 
e/(2 ms), where the mass of the baryon under investiga- 
tion appears, to the nuclear magneton, e/(2TOjv), where 



A. Charge radii 

To make contact with the extensive phenomenology of 
the field, our results for the electric form factors are ex- 
pressed in terms of charge radii. It is well known that the 
experimentally measured electric (and magnetic) form 
factor of the proton is described well by a dipole ansatz 
at small Q 2 



Qe{Q 2 ) 



Ge(0) 



(1 + Q 2 /to 2 ) 2 ' 



Q 2 >0. 



(5.1) 



where m characterises the size of the proton. This be- 
haviour has also been observed in recent lattice calcu- 
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TABLE II: Quark sector contributions to the electric form factors of the nucleon at Q 2 — 0.227(2) GeV 2 . Sector contributions 
are for a single quark having unit charge. The fit windows are selected using the criteria outlined in Sec. 1111131 



(GeV 2 ) 



fit value 



fit window 



X 2 /dof 



fit value 



dp 
fit window 



X 2 /dof 



0.9972(55) 


0.798(5) 


21-25 


1.02 


0.805(4) 


21-25 


2.36 


0.8947(54) 


0.789(5) 


21-25 


0.89 


0.796(5) 


21-25 


2.53 


0.7931(53) 


0.779(6) 


21-25 


0.64 


0.788(5) 


21-25 


2.21 


0.6910(35) 


0.768(6) 


21-25 


0.86 


0.780(5) 


21-25 


1.57 


0.5925(33) 


0.756(7) 


21-25 


0.80 


0.769(6) 


21-25 


1.42 


0.4854(31) 


0.740(9) 


21-25 


0.62 


0.755(9) 


21-25 


1.19 


0.3795(31) 


0.725(10) 


19-23 


1.23 


0.741(11) 


19-23 


0.78 


0.2839(33) 


0.708(12) 


19-23 


1.37 


0.723(14) 


19-23 


1.31 


0.2153(35) 


0.693(15) 


19-23 


0.82 


0.700(20) 


19-23 


1.23 


0.1384(43) 


0.682(17) 


16-20 


1.02 


0.678(25) 


16-20 


0.89 


0.0939(44) 


0.666(25) 


16- 19 


1.47 


0.644(38) 


16- 19 


1.28 



TABLE III: 
are for a sir 


Quark sector contributions to the electric form factors of Y. 
tgle quark having unit charge. The fit windows are selected 


', baryc 
using 


sirs at Q 2 = 0.227(2) 
the criteria outlined 


GeV 2 
in Sec. 


. Sector contributions 

linnl 


ml 

(GeV 2 ) 


fit value 


u E or d E 
fit window 




X 2 /dof 




fit value 




fit window 


X 2 /dof 



0.9972(55) 


0.793(6) 


21-25 


0.91 


0.759(6) 


21-25 


1.70 


0.8947(54) 


0.785(7) 


21-25 


0.86 


0.758(7) 


21-25 


1.71 


0.7931(53) 


0.776(7) 


21-25 


0.66 


0.757(8) 


21-25 


1.59 


0.6910(35) 


0.766(6) 


21-25 


1.00 


0.757(6) 


21-25 


1.40 


0.5925(33) 


0.755(7) 


21-25 


0.90 


0.756(7) 


21-25 


1.36 


0.4854(31) 


0.740(9) 


21-25 


0.62 


0.755(9) 


21-25 


1.19 


0.3795(31) 


0.726(10) 


19-23 


1.46 


0.754(10) 


19-23 


0.37 


0.2839(33) 


0.711(12) 


19-23 


1.78 


0.753(11) 


19-23 


0.58 


0.2153(35) 


0.700(14) 


19-23 


0.73 


0.752(13) 


19-23 


0.39 


0.1384(43) 


0.680(18) 


19-21 


0.73 


0.754(17) 


19-21 


0.18 


0.0939(44) 


0.670(23) 


19-23 


1.30 


0.750(26) 


19-21 


1.40 



TABLE IV: 
single quark 


Quark sector contributions to the electric form factors of A at Q 2 = 0.227(2) 
having unit charge. The fit windows are selected using the criteria outlined 


GeV 2 
in Sec. 


. Sector contributions 

utttjI 


are for a 


ml 
(GeV 2 ) 


fit value 


ma or d,A 
fit window 


X 2 /dof 




fit value 




SA 

fit window 


X 2 /dof 



0.9972(55) 


0.803(5) 


21-25 


1.20 


0.745(8) 


21-25 


0.64 


0.8947(54) 


0.794(6) 


21-25 


1.23 


0.744(9) 


21-25 


0.58 


0.7931(53) 


0.785(7) 


21-25 


1.06 


0.744(10) 


21-25 


0.54 


0.6910(35) 


0.775(6) 


21-25 


1.17 


0.738(8) 


21-25 


0.55 


0.5925(33) 


0.765(7) 


21-25 


1.12 


0.737(9) 


21-25 


0.48 


0.4854(31) 


0.750(8) 


21-25 


1.02 


0.735(10) 


21-25 


0.49 


0.3795(31) 


0.736(9) 


19-23 


0.94 


0.734(11) 


19-23 


0.88 


0.2839(33) 


0.720(11) 


19-23 


1.17 


0.730(12) 


19-23 


1.05 


0.2153(35) 


0.704(13) 


19-23 


1.23 


0.727(13) 


19-23 


0.55 


0.1384(43) 


0.694(13) 


16- 19 


3.79 


0.727(13) 


16- 17 


0.62 


0.0939(44) 


0.686(13) 


16- 17 


2.43 


0.729(14) 


16- 17 


0.29 
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TABLE V: Quark sector contributions to the electric form factors of 3 baryons at Q 2 = 0.227(2) GeV 2 . Sector contributions 
are for a single quark having unit charge. The fit windows are selected using the criteria outlined in Sec. 1111 HI 



(GeV 2 



fit value 



fit window 



X 2 /dof 



fit value 



Us or ds 
fit window 



X 2 /dof 



0.9972(55) 


0.747(9) 


21-25 


0.34 


0.804(8) 


21-25 


1.60 


0.8947(54) 


0.747(9) 


21-25 


0.36 


0.794(8) 


21-25 


1.53 


0.7931(53) 


0.746(10) 


21-25 


0.37 


0.785(9) 


21-25 


1.52 


0.6910(35) 


0.742(8) 


21-25 


0.54 


0.778(7) 


21-25 


1.38 


0.5925(33) 


0.741(8) 


21-25 


0.55 


0.768(8) 


21-25 


1.24 


0.4854(31) 


0.740(9) 


21-25 


0.62 


0.755(9) 


21-25 


1.19 


0.3795(31) 


0.739(9) 


19-23 


0.70 


0.740(10) 


21-25 


1.45 


0.2839(33) 


0.738(10) 


19-23 


1.18 


0.723(13) 


21-25 


1.22 


0.2153(35) 


0.736(10) 


19-23 


1.42 


0.709(16) 


21-25 


0.81 


0.1384(43) 


0.733(10) 


19-23 


0.52 


0.690(19) 


20-23 


0.71 


0.0939(44) 


0.725(11) 


19-23 


1.21 


0.672(22) 


20-23 


0.59 



TABLE VI: 
are for a sin 


Quark sector contributions to the magnetic form factors of the nucleon at Q 2 = 0.227(2) GeV 2 
gle quark having unit charge. The fit windows are selected using the criteria outlined in Sec. 


. Sector contributions 

IiTTTjI 


m 2 (GeV 2 ) 


fit value 


u p (un) 
fit window 


X 2 /dof 


fit value 


dp (un) 
fit window 


X 2 /dof 



0.9972(55) 


0.765(12) 


19-23 


1.78 


-0.295(7) 


18-22 


0.61 


0.8947(54) 


0.785(14) 


19-23 


1.33 


-0.298(8) 


18-22 


0.51 


0.7931(53) 


0.804(16) 


19-23 


1.01 


-0.301(9) 


18-22 


0.43 


0.6931(51) 


0.817(13) 


19-23 


1.00 


-0.301(8) 


18-22 


0.91 


0.5944(51) 


0.838(15) 


19-23 


0.73 


-0.304(10) 


18-22 


0.79 


0.4869(50) 


0.861(20) 


19-23 


0.64 


-0.306(12) 


18-22 


0.86 


0.3795(31) 


0.893(24) 


17-21 


0.14 


-0.314(14) 


16- 18 


1.64 


0.2839(33) 


0.932(31) 


17-21 


0.14 


-0.313(19) 


16-19 


1.24 


0.2153(35) 


0.967(42) 


17-21 


0.61 


-0.313(31) 


16-20 


0.53 


0.1384(43) 


1.034(52) 


16-20 


1.12 


-0.309(40) 


15-19 


0.49 


0.0939(44) 


1.024(72) 


15- 17 


0.82 


-0.335(54) 


15-17 


2.44 



TABLE VII: Quark sector contributions to the magnetic form factors of E baryons at Q — 0.227(2) GeV . Sector contributions 
are for a single quark having unit charge. The fit windows are selected using the criteria outlined in Sec. II1IDI 



m 2 (GeV 2 



fit value 



u s or d s (/j, n ) 
fit window 



xVdof 



fit value 



ss (un) 
fit window 



X 2 /dof 



0.9972(55) 


0.775(14) 


19-23 


1.22 


-0.316(10) 


18-22 


0.73 


0.8947(54) 


0.793(16) 


19-23 


1.10 


-0.314(11) 


18-22 


0.63 


0.7931(53) 


0.810(18) 


19-23 


0.96 


-0.312(11) 


18-22 


0.55 


0.6910(35) 


0.821(14) 


19-23 


0.98 


-0.309(9) 


18-22 


0.94 


0.5925(33) 


0.840(16) 


19-23 


0.82 


-0.308(10) 


18-22 


0.86 


0.4854(31) 


0.861(20) 


19-23 


0.64 


-0.306(12) 


18-22 


0.86 


0.3795(31) 


0.886(23) 


17-21 


0.19 


-0.308(13) 


16-18 


0.42 


0.2839(33) 


0.914(27) 


17-21 


0.26 


-0.310(15) 


16-18 


0.03 


0.2153(35) 


0.941(32) 


17-21 


0.61 


-0.317(19) 


16-18 


0.04 


0.1384(43) 


0.964(33) 


15-20 


1.57 


-0.317(24) 


16-18 


0.92 


0.0939(44) 


0.969(41) 


15- 17 


0.11 


-0.322(28) 


16- 18 


1.08 
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TABLE VIII: Quark sector contributions to the magnetic form factors of A at Q 2 = 0.227(2) GeV 2 . Sector contributions are 
for a single quark having unit charge. The fit windows are selected using the criteria outlined in Sec. lIIIDl 



(GeV 2 



fit value 



ua or dA (mjv) 
fit window 



X 2 /dof 



fit value 



SA ((J,n) 

fit window 



X 2 /dof 



0.9972(55) 


0.069(9) 


19-23 


1.26 


1.200(22) 


18-22 


1.38 


0.8947(54) 


0.072(9) 


19-23 


1.20 


1.210(22) 


18-22 


1.17 


0.7931(53) 


0.075(10) 


19-23 


1.07 


1.214(24) 


18-22 


0.99 


0.6910(35) 


0.079(8) 


19-23 


1.08 


1.217(17) 


18-22 


0.95 


0.5925(33) 


0.083(9) 


19-23 


1.00 


1.228(18) 


18-22 


0.80 


0.4854(31) 


0.087(11) 


19-23 


0.88 


1.241(21) 


18-22 


0.75 


0.3795(31) 


0.091(12) 


17-21 


0.55 


1.259(22) 


16-20 


0.59 


0.2839(33) 


0.097(13) 


17-21 


0.29 


1.277(25) 


16-20 


0.42 


0.2153(35) 


0.099(17) 


17-21 


0.36 


1.293(28) 


16-20 


0.42 


0.1384(43) 


0.105(18) 


15- 16 


1.79 


1.308(33) 


16- 18 


0.94 


0.0939(44) 


0.105(22) 


15- 16 


0.19 


1.315(37) 


15- 16 


1.98 



TABLE IX: Ql 
are for a single 


ark sector contributions to the magnetic form factors of 
quark having unit charge. The fit windows are selected 


3 bary 
using 


ons at Q 2 = 0.227(2) 
the criteria outlined 


GeV 2 . Sector contributions 
in Sec.liTTTjl 


m 2 (GeV 2 ) 


fit value 


ss Qltjv) 
fit window 




X 2 /dof 




fit value 


Us 


or cfe (A* at) 
fit window 


X 2 /dof 



0.9972(55) 


0.846(27) 


19-23 


0.60 


-0.290(10) 


18-22 


0.19 


0.8947(54) 


0.849(27) 


19-23 


0.67 


-0.294(11) 


18-22 


0.21 


0.7931(53) 


0.851(28) 


19-23 


0.75 


-0.299(12) 


18-22 


0.21 


0.6931(51) 


0.855(19) 


19-23 


0.41 


-0.300(10) 


18-22 


0.59 


0.5944(51) 


0.858(19) 


19-23 


0.49 


-0.303(11) 


18-22 


0.68 


0.4869(50) 


0.861(20) 


19-23 


0.64 


-0.306(12) 


18-22 


0.86 


0.3795(31) 


0.866(21) 


16-20 


0.84 


-0.311(13) 


16-20 


1.77 


0.2839(33) 


0.871(21) 


16-20 


0.78 


-0.316(14) 


16-20 


1.66 


0.2153(35) 


0.875(22) 


16-20 


1.04 


-0.322(15) 


16-20 


0.77 


0.1384(43) 


0.879(22) 


16-20 


1.16 


-0.328(17) 


15-19 


0.86 


0.0939(44) 


0.879(23) 


16-20 


0.55 


-0.333(18) 


15-19 


1.00 



lations [12} where many momentum transfers have been 
considered. Using this observation, together with 



( 



i> 



-6 



w Ge(q2 



(5.2) 



? 2 =o 



we arrive at an expression which allows us to calcu- 
late the electric charge radius of a baryon using our 
two available values of the Sach's electric form factor 
(MQmiJ. M0)), namely 
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(r|) __ 
&s (0) ~ Q 2 



I Ge(0) 
Ge(Q 2 ) 



(5.3) 



While Eq. I|5.1|l is suitable for a charged baryon, al- 
ternative forms must be considered for neutral baryons 
where Ge(0) = 0. 



However, we have direct access to the charge distribu- 
tions of the individual quark sectors, a subject receiving 
tremendous experimental attention in the search for the 
role of hidden flavour in baryon structure. In this case 
Eq. i|5.3|) may be applied to each quark sector providing 
an opportunity to determine the charge radii on a sector 
by sector basis. 

For neutral baryons it becomes a simple matter to con- 
struct the charge radii by first calculating the charge 
radii for each quark sector. These quark sectors are then 
combined using the appropriate charge and quark num- 
ber factors as described in Sec. II V Hi to obtain the total 
baryon charge radii. Indeed, all baryon charge radii, in- 
cluding the charged states, are calculated in this manner. 

Tables [X] to IXIIII provide the electric charge radii of 
the octet baryons and their quark-sector contributions 
normalised to the case of single quarks with unit charge. 
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TABLE X: Nucleon electric charge radii squared. Quark-sector contributions are indicated for single quarks having unit charge. 
Baryon charge states are also summarised. Values for m^ and (r 2 ) are in units of GeV 2 and fm 2 respectively. 



(In 



P 



0.9972(55) 
0.8947(54) 
0.7931(53) 
0.6910(35) 
0.5925(33) 
0.4854(31) 
0.3795(31) 
0.2839(33) 
0.2153(35) 
0.1384(43) 
0.0939(44) 



0.243(7) 
0.256(8) 
0.270(9) 
0.288(9) 
0.307(11) 
0.332(14) 
0.358(17) 
0.389(22) 
0.416(27) 
0.437(31) 
0.467(48) 



0.231(6) 
0.245(7) 
0.257(8) 
0.270(8) 
0.286(10) 
0.309(14) 
0.332(17) 
0.363(24) 
0.403(36) 
0.445(46) 
0.510(77) 



0.247(8) 
0.259(9) 
0.273(10) 
0.294(10) 
0.314(12) 
0.340(16) 
0.367(19) 
0.397(24) 
0.420(29) 
0.435(32) 
0.452(53) 



-0.007(3) 
-0.007(3) 
-0.008(4) 
-0.012(4) 
-0.014(5) 
-0.015(7) 
-0.017(9) 
-0.017(12) 
-0.008(16) 
0.005(21) 
0.029(41) 



TABLE XI: E electric charge radii squared. Quark-sector contributions are indicated for single quarks having unit charge. 
Baryon charge states are also summarised where absolute values of the E~ results are reported. Values for m 2 and (r 2 ) are in 
units of GeV and fm respectively. 



ml 


we or ds 


SE 


E+ 


E° 


E~ 


0.9972(55) 


0.249(9) 


0.301(10) 


0.232(10) 


-0.017(3) 


0.266(9) 


0.8947(54) 


0.261(10) 


0.302(11) 


0.249(12) 


-0.013(3) 


0.275(9) 


0.7931(53) 


0.276(11) 


0.304(12) 


0.266(15) 


-0.010(3) 


0.285(11) 


0.6910(35) 


0.291(10) 


0.304(10) 


0.286(11) 


-0.005(2) 


0.295(9) 


0.5925(33) 


0.309(11) 


0.306(11) 


0.310(13) 


0.001(3) 


0.308(11) 


0.4854(31) 


0.332(14) 


0.309(14) 


0.340(16) 


0.008(4) 


0.324(13) 


0.3795(31) 


0.356(16) 


0.310(15) 


0.371(19) 


0.015(4) 


0.341(15) 


0.2839(33) 


0.382(20) 


0.312(18) 


0.405(23) 


0.023(5) 


0.359(18) 


0.2153(35) 


0.401(25) 


0.315(21) 


0.430(29) 


0.029(6) 


0.372(22) 


0.1384(43) 


0.438(33) 


0.311(26) 


0.480(38) 


0.042(8) 


0.395(28) 


0.0939(44) 


0.456(43) 


0.318(41) 


0.503(54) 


0.046(15) 


0.410(37) 



1. Quenched chiral perturbation theory 

The effective field theory formalism of quenched chiral 
perturbation theory (QxPT) predicts significant contri- 
butions to the charge radii which have their origin in 
virtual meson-baryon loop transitions. These loops give 
rise to contributions which have a non-analytic depen- 
dence on the quark mass or squared pion mass. While 
the absence of sea-quark loops generally acts to suppress 
the magnitude of the coefficients of these terms (and oc- 
casionally the sign is reversed), there are several channels 
in which these contributions remain significant. 

The leading non-analytic (LNA) and next-to-leading 
non-analytic (NLNA) behaviour of charge distribution 
radii in full QCD are 



v4)" J^2j2 E 



5p i log[^)~10 r 3' i g(m i ,A,n) 



-c + c 2 mf + C4 m i . 



(5.4) 



Here the sum over i includes the it and K pseudoscalar 
mesons. The contributions of the various charge states of 
these mesons are contained in the coefficients f3 and /3', 
reflecting electric charge and SU(3) axial couplings, D, 
F and C. In quenching the theory, the coefficients j3 and 
j3' are modified to reflect the absence of sea-quark loops. 

The first term arises from octet baryon to octet-baryon 
- meson transitions. Thus charge radii are charac- 
terised by a logarithmic divergence |35j in the chiral limit 
{ml — > 0). In this simple form, the mass splittings be- 
tween baryon octet members is neglected. 

The second term of Eq. I|5.4|) arises from octet baryon 
to decuplet-baryon - meson transitions. As the splitting 
between the baryon octet and decuplet does not vanish 
in the chiral limit, the mass splitting, A = Ma — Mn, be- 
tween the nucleon and A for example, plays an important 
role. The function Q(mi,A,,a) is 



S(m,A, M )=log(-i- 
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TABLE XII: A electric charge radii squared. Quark-sector contributions are indicated for single quarks having unit charge. 
The baryon charge state is also provided. Values for m 2 and (r 2 ) are in units of GeV 2 and fm 2 respectively. 



ma or di\ 



»\ 



A" 



0.9972(55) 
0.8947(54) 
0.7931(53) 
0.6910(35) 
0.5925(33) 
0.4854(31) 
0.3795(31) 
0.2839(33) 
0.2153(35) 
0.1384(43) 
0.0939(44) 



0.235(8) 
0.248(9) 
0.262(10) 
0.276(9) 
0.293(10) 
0.316(13) 
0.340(15) 
0.367(18) 
0.395(23) 
0.414(24) 
0.428(25) 



0.322(13) 
0.323(14) 
0.324(15) 
0.335(13) 
0.336(14) 
0.340(16) 
0.343(18) 
0.350(19) 
0.357(22) 
0.356(22) 
0.354(23) 



-0.029(3) 

-0.025(3) 

-0.021(3) 

-0.020(3) 

-0.014(3) 

-0.008(4) 

-0.001(4) 

0.006(5) 

0.013(6) 

0.019(6) 

0.025(7) 



TABLE XIII: — electric charge radii squared. Quark-sector contributions are indicated for single quarks having unit charge. 



Baryon charge states are also summarised where absolute values of the ! 
units of GeV and fm respectively. 



results are reported. Values for m^ and (r ) are in 



us or ds 



0.9972(55) 
0.8947(54) 
0.7931(53) 
0.6910(35) 
0.5925(33) 
0.4854(31) 
0.3795(31) 
0.2839(33) 
0.2153(35) 
0.1384(43) 
0.0939(44) 



0.319(14) 
0.320(15) 
0.322(16) 
0.329(13) 
0.330(14) 
0.332(14) 
0.334(15) 
0.337(15) 
0.340(16) 
0.345(16) 
0.358(18) 



0.235(11) 
0.249(12) 
0.262(13) 
0.273(11) 
0.288(12) 
0.309(14) 
0.333(17) 
0.361(21) 
0.385(27) 
0.419(34) 
0.451(41) 



0.290(12) 
0.296(13) 
0.301(14) 
0.310(11) 
0.316(12) 
0.324(13) 
0.334(14) 
0.345(16) 
0.355(18) 
0.370(19) 
0.389(23) 



-0.056(7) 
-0.047(7) 
-0.039(8) 
-0.037(6) 
-0.028(7) 
-0.015(7) 
-0.0003(80) 
0.016(10) 
0.029(13) 
0.049(19) 
0.062(22) 



V& 2 



:10g 



A-v/aX 



mt 



A + VA 1 



(5-5) 



As the tadpole graph contributing to the LNA term 
of charge radii in full QCD vanishes in quenched QCD 
[2a|, the coefficients /3 and /?' for charge radii are iden- 
tical to those for magnetic moments in quenched QCD 
[14], [HI- Figure [7| displays the non- analytic contributions 
from Q%PT as given in Eq. I|5.4[l . plotted for the sample 
case of the proton. In this case, the values of (3 and (3' 
are — |D 2 and — \C 2 respectively [lj,[lla]. Here the axial 
couplings D and C are related by C = 2D and D is taken 
as 0.76. The scale /i 2 is taken to be 1 GeV 2 and serves 
only to define cq. 

Because these non-analytic contributions are comple- 
mented by terms analytic in the quark mass or pion-mass 
squared, the slope and curvature at large m 2 of these 
contributions is not significant. What is significant is the 
curvature at small m 2 and we see that this curvature 
is dominated by the LNA term. Here there is no mass 
splitting to mask the effects of dynamical chiral symme- 



try breaking. Thus, we will examine the extent to which 
our simulation results are consistent with the LNA be- 
haviour of QxPT. 

The coefficient (3 is related to the coefficient of the 
leading non-analytic (LNA) contribution to the magnetic 
moment, \i y i a the relation [l4| 







m N 

MS 



X 



(5.6) 



The coefficients \ have been determined for octet 
baryons and their individual quark sector contributions 
in Ref. J14J and numerical values are reproduced in Ta- 
bles ]XIV\ and |^V] for ready reference. 

Since m 2 < 1 GeV 2 in our simulations, the logarith- 
mic term is negative for all quark masses considered here. 
Hence, the charge radius will exhibit a logarithmic diver- 
gence in chiral limit to either positive or negative infinity, 
depending on the whether (3 (or x) is negative or positive 
respectively. 

In the quenched approximation, the flavour-singlet rf 
meson remains degenerate with the pion and makes im- 
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FIG. 7: The leading (upper curve) and next-to-leading (lower 
curve) non-analytic contributions to the charge radius of the 
proton as given by quenched chiral perturbation theory in 
Eq. |JOJ. 



portant contributions to quenched chiral non-analytic be- 
haviour. The neutrality of its charge prevents it from 
contributing to the coefficients of Tables IXIVI and IXVI 
However, the double hair-pin diagram in which the vec- 
tor current couples to the virtual baryon intermediary 
does give rise to chirally-singular behaviour. However 
the relatively small couplings render these contributions 
small at the quark masses probed here. 



2. Quark sector charge radii 

We begin with an examination of the quark contribu- 
tions to baryon charge radii. The results are reported for 
single quarks of unit charge. Of particular interest are 
the contributions of similar quarks experiencing differ- 
ent environments. Traditionally, quark models of hadron 
structure neglected such environment sensitivity. How- 
ever, such environment sensitivity is manifest in chiral 
effective field theory. The finite kaon mass in the chi- 
ral limit renders the kaon's contributions to curvature 
almost trivial relative to the pion. 

Figure[S]displays the charge radii of the w-quark distri- 
bution in the proton and compares this with the u-quark 
distribution in E + . The 5'[/(3)-flavour limit is manifest 
at m 2 . ~ 0.5 GeV 2 . The replacement of a d-quark in the 
proton, by an s quark in E + gives rise to only a small 
environment sensitivity in the u-quark properties. 

Referring to the chiral coefficients of Table IXIVI the 
negative value of x for u p indicates that the charge ra- 
dius of the u quark distribution in the proton should 
diverge to positive infinity in the chiral limit. A physical 
understanding of this is made obvious by considering the 
virtual transition p — > mr + , which at the quark level can 
be understood as (uud) — > (udd)(du). In the chiral limit, 
the 7r + carries the u quark to infinity such that u-quark 



TABLE XIV: Coefficients, x, providing the LNA contribution 
to baryon magnetic moments and charge radii in quenched 
QCD. Coefficients for magnetic moments in full QCD are 
also indicated. Here the coefficients for quark sector contri- 
butions to baryon properties are indicated for quarks having 
unit charge. Note that u p for example denotes the coefficient 
for the two u quarks of the proton, each of which have unit 
charge. Intermediate (Int.) meson-baryon channels are indi- 
cated to allow for S'C/(3)-flavour breaking in both the meson 
and baryon masses. The coefficients are calculated from the 
expressions of Ref. p3 with the axial couplings F — 0.50 and 
D = 0.76 with f n = 93 MeV. 



tip | &n 



dp | tin 



u s + I d s - 



d E+ | u E - 



-ss 



u^o | d E o 



"A I d\ 



•SA 



m s o | d s - 



d=n | fe- 



int. 



Ntt 
AK 
T,K 

Ntt 
T.K 

AK 

T,K 



A-k 
NK 
SK 

Ett 
A-k 
NK 

NK 
5K 
Etj s 

Stt 
Att 
NK 
EK 



Y.TT 

A-qi 
NK 
5K 

Arj s 
NK 
EK 



■=,7T 

AK 
EK 
Q.K 

En 
AK 
EK 

AK 
T,K 
Q.K 

Erj s 



Full QCD 



-6.87 
-3.68 
-0.15 

+6.87 
-0.29 

+3.68 
+0.44 



-2.16 
-1.67 


-6.87 

+2.16 
+1.67 
+0.29 

-0.29 

+6.87 






+0.15 
-3.43 






+3.68 
-0.40 


-7.36 
+0.79 



-0.29 



+6.87 



+0.29 
+0.40 
+3.43 

-0.40 
-10.3 







Quenched QCD 



-3.33 




+3.33 










-0.29 
-3.04 





+0.29 

+3.04 






-0.15 
-1.52 






+1.23 
+0.44 



-2.45 
-0.88 




-0.40 
+3.43 
+0.29 






+0.40 

-3.43 

-0.29 





Hi 



TABLE XV: Coefficients, \, providing the LNA contribution 
to baryon magnetic moments and charge radii in quenched 
QCD. Coefficients for magnetic moments in full QCD are also 
indicated. Intermediate (Int.) meson-baryon channels are in- 
dicated to allow for S'(7(3)-flavour breaking in both the meson 
and baryon masses. The coefficients are calculated from the 
expressions of Ref. [l4J | with the axial couplings F — 0.50 and 
D = 0.76 with /„ = 93 MeV. 



Baryon 



Channel 



Full QCD 



A 



Nn 
AK 
SIC 

Ayr 
AK 
EK 

Ett 

A7T 

NK 
EK 

T, Vs 

Ett 
An 
NK 
EK 
Et? s 

Ett 
Aty 
NK 
EK 
Et? s 

Ett 
Arn 
NK 
EK 
Arj s 

En 
AK 
T,K 

file 

Er/s 

En 
AK 
T,K 
QK 

Erj s 



-6.87 


-3.68 


-0.15 


+6.87 





-0.29 


-2.16 


-1.67 





-6.87 











+0.15 


-3.43 





+2.16 


+1.67 


+0.29 














+3.68 


-0.40 





-0.29 





+6.87 








+0.29 


+0.40 


+3.43 









Quenched QCD 



-3.33 








+3.33 














-0.29 


-3.04 











-0.15 


-1.52 


























+1.23 


+0.44 








-0.40 


+3.43 


+0.29 
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FIG. 8: Electric charge radii of the u-quark distribution in the 
proton, Up, and E + , u E +, as a function of m^. representing 
the quark masses considered in the simulation. The results 
for Up are offset for clarity. 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.B 0.9 1.0 1.1 
m w z (GeV 2 ) 

FIG. 9: Charge radii of u n and u s o as a function of m\. The 
data for the u n are offset for clarity. 



charge distribution radius in the proton diverges. 

In the case of the u quark in S + , the coefficient of the 
logarithmic divergence is zero in the 7r channel and hence 
no divergence is expected. While there is a significant co- 
efficient for transitions to 'EK, the increased mass of the 
S baryon makes this channel unfavourably suppressed. 

The results in Fig.[S]for u p exhibit an upward trend and 
increasing curvature with reducing quark mass. The m s + 
rises more slowly. Hence these results are in qualitative 
agreement with the LNA expectations of chiral effective 
field theory. 

FigureEldisplays the electric charge radii of the u quark 
in the neutron (u n ) and in S° (uh°) as a function of mj. 
Here we observe that the charge radii of u n and u s o are 
nearly equal at heavy quark masses, but in the chiral 
limit they differ. The light d-quark environment of the u 
quark in the neutron provides enhanced chiral curvature 
as the chiral limit is approached. 
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However, the true nature of the underlying physics is 
much more subtle. From Table 1X1 VI we see that the 
quenched coefficient (last column) for the u quark in the 
neutron is positive in the 7r channel, from which we can 
deduce that the charge radius should actually diverge to 
negative infinity in the chiral limit. 

Physically this can be understood by looking at the 
quark contributions to the virtual transition n — > pTt~ 
which gives rise to this divergence. In this case one has 
(ddu) — ► (duu)(ud). In the chiral limit, the mass of the 
pion approaches zero such that the tt~ carries a u quark 
to infinity. Since the d quark is ignored while calculating 
the u quark contribution {i.e. the electric charge of the 
d quark may be thought of as zero) , the entire charge of 
the pion comes from the u quark, thus taking the u-quark 
charge distribution radius to negative infinity. However, 
Fig. El shows no such trend. 

The coefficient \ f° r u s° is zero in the 7r channel, indi- 
cating that there should be no logarithmic divergence in 
the chiral limit. However it does have a substantial posi- 
tive coefficient in the favourable T,K channel, indicating 
the possibility of downward curvature as the chiral limit 
is approached. Again, Fig. shows no hint of downward 
curvature. 

While the statistical error bars are sufficiently large 
to hide such a turn over, there are other interpretations. 
One possibility is that we are not yet in the true chiral 
regime where such physics is manifest. Indeed, the di- 
vergence of the u-quark charge distribution to negative 
infinity may only reveal itself at quark masses lighter than 
the physical quark masses. 

Alternatively, one might regard this particular case 
to be somewhat exceptional. It is the only channel in 
which chiral-loop physics is expected to oppose the nat- 
ural broadening of a distribution's Compton wavelength. 
On the lattice, the finite volume restricts the low mo- 
menta of the effective field theory to discrete values. It 
may be that this lattice artifact prevents one from build- 
ing up sufficient strength in the loop integral to counter 
the Compton broadening. In this case it would be im- 
possible to observe the divergence of u n — » — oo at any 
quark mass. It will be interesting to resolve this discrep- 
ancy with quantitative effective field theory calculations. 

Figure [IJ]| reports our results for the charge distribu- 
tion radius of a u quark in A as a function of mj. The 
chiral coefficient for this is zero in the it channel and 
hence no divergence is expected. However, there is sig- 
nificant strength for downward curvature in the energet- 
ically favourable NK channel. Indeed, the approach to 
the chiral limit is remarkably linear and contrasts the up- 
ward curvature observed for other light quark flavours. 
Hence our results are in qualitative agreement with the 
expectations of QxPT. 

Figure ITT1 illustrates the charge distribution radius of 
strange quarks in A , S° and S°. In our simulations the 
strange quark mass is held fixed and therefore any vari- 
ation observed in the results is purely environmental in 
origin. All three distributions suggest a gentle depen- 
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dence on the mass of the environmental light quarks. 

However, the environmental flavour-symmetry depen- 
dence of the strange quark distributions is absolutely 
remarkable. When the environmental quarks are in an 
isospin state in the A, the strange quark distribution 
is broad. On the other hand, when the environmental 
quarks are in an isospin 1 state in X baryons, the distri- 
bution radius is significantly smaller. 

In the case of strange quark distributions, the LNA 
contributions are exclusively from transitions involving 
the kaon. Therefore significant curvature is not expected. 
On the other hand, one might expect broader distribu- 
tions in cases where a virtual transition is possible in 



quenched QCD. Referring to Table IXIVI one sees that 
both sa and Ss have strong transitions to the energeti- 
cally favourable KN and KE channels respectively. The 
coefficients are negative such that the virtual transitions 
will act to enhance the charge distributions. This is not 
the case for ss where the sign is positive and the tran- 
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sition is to the energetically unfavoured KE channel. In 
summary, QxPT suggests the charge distributions for s\ 
and sh will be larger than for S£. This is exactly what 
is observed in Figure ITT1 



3. Baryon charge radii 

The flavour-symmetry dependence of uds-quark distri- 
butions in E° and A is particularly manifest in Fig. ^| 
Here the interplay between the light-quark sector with 
effective charge 1/3 and the strange sector with charge 
— 1/3 is revealed. 

At the SU{2>) flavour-symmetric point (m^ ~ 
0.5 GeV 2 ) where the strange and light quarks have the 
same mass and the A and £ are degenerate in mass, nei- 
ther charge radius is zero. This very nicely reveals differ- 
ent charge distributions for the quark sectors described 
in the previous section. 

In constituent quark models, this flavour dependence 
would be described in terms of spin-dependent forces. 
In the A where a scalar diquark can form between the 
non-strange pair, the charge radius is dominated by the 
broader strange-quark distribution at the S't/(3)-flavour 
symmetric point. This is contrasted by the E° where 
scalar-diquark pairing would occur between strange and 
non-strange quarks, acting to constrict the strange quark 
distribution in S as seen in Fig. ^] In addition, hyper- 
fine repulsion in the non-strange quark sector leads to 
a broader distribution for the light quark sector as indi- 
cated in Tables OTA and IXTTl As compelling as this discus- 
sion is, this line of reasoning suggests the decuplet baryon 
states should have broader quark distributions [20J as 
scalar-diquark clusters do not dominate there. However, 
preliminary results from an analysis of decuplet baryon 
structure on the same lattice configurations explored here 
|37j . do not reveal broader quark distributions. For this 
reason, we consider our discussion of virtual transitions 
in the context of effective field theory in the previous sec- 
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tion to be a more relevant description of the underlying 
physics. 

Ultimately, as the chiral limit is approached, the light 
quark distribution broadens and dominates the charge 
radii for both baryons. However, the charge distribution 
of the £ is much broader and reflects our discussion of 
the quark sector contributions. In particular, the LNA 
contributions of Qx?T act to suppress the distribution 
of ma and enhance s\, whereas the LNA contributions 
to E° are relatively suppressed either by having small 
coefficients or having energetically unfavourable transi- 
tions in the kaon channel. This suppression of ma and 
enhancement sa combines to give a strong net effect of 
suppressing the charge radius of the A . 

The hyperon charge states, S~ and S~ have chiral co- 
efficients which vanish in quenched QCD. Similarly, E + 
has no contributions from virtual pion transitions. The 
one case, where there is a substantial coefficient, is sup- 
pressed energetically. Figure IT31 displays our simulation 
results for the electric charge distribution radii of these 
hyperons as a function of to 2 . The ordering of the charge 
radii as the chiral limit is approached is explained by the 
more localized strange quark distribution. 

Figure [21 compares the charge radii of E + with the 
proton. The charge radii of these baryons match at the 
SU(3) flavour limit where m 2 ~ 0.5 GeV 2 as required. 
As the chiral limit is approached, the smaller charge dis- 
tribution of the heavier negatively-charged strange quark 
acts to make the £+ larger. This is manifest in the sim- 
ulation results. 

While the £+ is not expected to display chiral curva- 
ture, the proton charge radius presents one of the more 
favourable opportunities to observe a hint of the loga- 
rithmic divergence to be encountered in the chiral and 
infinite-volume limits of quenched QCD. However, there 
is no hint of chiral curvature in favour of the proton over 
the S+. 

The origin of this discrepancy is once again traced to 
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the singly-represented u-quark in the neutron, or more 
specifically in this case, the singly-represented d-quark 
in the proton. As highlighted in the discussion surround- 
ing Fig. [SI there is a hint of increased curvature for the 
doubly-represented w-quark in the proton over that in S, 
in accord with chiral effective field theory. But this is hid- 
den in the proton charge radius due to the absence of the 
anticipated curvature of the singly-represented quark in 
the nucleon, as highlighted in the discussion surrounding 
Fig.H 

Similarly, the ultimate divergence of the neutron 
charge radius to negative infinity via n — > pir~ is not 
yet manifest. Rather a crossing of the central values 
into positive values of squared charge radii is revealed 
in Fig. ED Still, the statistical errors remain consistent 
with negative values. 

The crossing of the central values of the squared neu- 
tron charge radius into positive values led us to further 
examine our selection of fit regime in our correlation func- 
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FIG. 16: The proton charge radius is compared with previous 
state of the art lattice simulation results in quenched QCD. 
The solid squares indicate current lattice QCD results with 
FLIC fermions. The stars indicate the lattice results of Q 
while the crosses indicate the results of [EJ , both of which use 
the standard Wilson actions for the gauge and fermion fields. 

tion analysis. Our concern is that noise in the correlation 
function may be distorting the fit. Hence, we have also 
considered fits including t = 15, immediately following 
the point-split current insertion centred at t = 14. While 
we prefer to allow some Euclidean time evolution follow- 
ing the current insertion, this systematic uncertainty is 
reflected in the asymmetric error bar of Fig. ED for the 
lightest two neutron charge radii. 

To summarise, we have explored the electric form fac- 
tors of the baryon octet and their quark sector contribu- 
tions at light quark masses approaching the chiral regime. 
The unprecedented nature of our quark masses is illus- 
trated in Fig. ED which compares the present results for 
the proton charge radius with the previous state of the 
art 0, 13 ■ Here the static quark potential has been used 
to uniformly set the scale among all the results. The 
small values of the early results are most likely due to 
the small physical lattice volumes necessitated at that 
time. The precision afforded by 400 20 3 x 40 lattices is 
manifest. 

We have discovered that all baryons having non- 
vanishing energetically-favourable couplings to virtual 
meson-baryon transitions tend to be broader than those 
which do not. This qualitative realisation provides a sim- 
ple explanation for the patterns revealed in our quenched- 
QCD simulations. 

Still, evidence of chiral curvature on our large-volume 
lattice is rather subtle in general and absent in the ex- 
ceptional case of the singly-represented quark in the neu- 
tron or S. In this case, it is thought that the restriction 
of momenta to discrete values on the finite-volume lat- 
tice prevents the build up of strength in the loop integral 
of effective field theory. Without sufficient strength, the 
Compton broadening of the distribution will not be coun- 
tered as the chiral limit is approached. 
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B. Magnetic moments 

The magnetic moment /i is provided by the magnetic 
form factor at Q 2 = 0, Gm(Q), with units of the natural 
magneton, /j,g = e/(2Mg), where Mg is the mass of the 
baryon 



M = Sm(0) 



2M b 



(5.7) 



While we could present a detailed discussion of the mag- 
netic form factors summarised in Sec. IIVBI a more in- 
teresting discussion of the results is facilitated via the 
magnetic moment where chiral non-analytic behaviour 
takes on a simple functional form and a vast collection of 
phenomenology is available to provide a context for our 
results. 

Since the magnetic form factors must be calculated at 
a finite value of momentum transfer, Q 2 , the magnetic 
moment must be inferred from our results, Gm(Q 2 ), ob- 
tained at the minimum non-vanishing momentum trans- 
fer available on our periodic lattice. The Q 2 dependence 
of lattice results from the QCDSF collaboration [12J are 
described well by a dipole. Phenomenologically this is a 
well established fact for the nucleon at low momentum 
transfers. 

However, we will take an even weaker approximation 
and assume only that the Q 2 dependence of the electric 
and magnetic form factors is similar, without stating an 
explicit functional form for the Q 2 dependence. This 
too is supported by experiment where the proton ratio 

Q (o*\ ~ •"■ f° r va l ues °f Q 2 similar to that probed here. 
In this case 



9m(P) 



Qm{Q 2 ) 
Ge(Q 2 ) 



MO). 



(5.8) 



The strange and light sectors of hyperons will scale dif- 
ferently, and therefore we apply Eq. (|5.8|) to the indi- 
vidual quark sectors. Octet baryon properties are then 
reconstructed as described in the discussion surround- 
ing Eq. (|4.1|l in Sec. IIVBI Results for baryon magnetic 
moments and their quark sector contributions are sum- 
marised in Tables IXVll through IXlXl 
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FIG. 17: Magnetic moment contributions of the u-quark sec- 
tor to the proton, u p , and the E + , u E +. The contributions 
Up are shifted right for clarity. 

is finite in the chiral limit, the rate of change of this con- 
tribution does indeed diverge in the chiral limit. The less 
singular nature of this contribution should allow its con- 
tributions to be observed at larger pion masses, making 
magnetic moments an excellent observable to consider 
in searching for evidence of chiral curvature. Kaon con- 
tributions take on the same form in the limit in which 
baryon mass splittings are neglected. 

As for the charge radii, negative values of x provide 
curvature towards more positive values as the chiral limit 
is approached, and vice versa for positive values of %. 

As emphasised earlier in our discussion of charge radii, 
the flavour-singlet rj' meson remains degenerate with the 
pion in the quenched approximation and makes impor- 
tant contributions to quenched chiral non-analytic be- 
haviour. The neutrality of its charge prevents it from 
contributing to the coefficients of Tables IXIVI and IXVI 
However, the double hair-pin diagram in which the vec- 
tor current couples to the virtual baryon intermediary 
does give rise to a logarithmic divergence in baryon mag- 
netic moments. However the relatively small couplings of 
the rj render these contributions negligible at the quark 
masses probed here |3q . 



1. Quenched chiral perturbation theory 

As for the charge radii, it is interesting to compare our 
results with the LNA and NLNA terms of xPT which 
survive to some extent in QxPT. As for the char ge r adii, 
the NLNA contributions provide little curvature (38| and 
we turn our attention to the LNA contributions |14j . 
These LNA contributions to baryon magnetic moments 
have their origin in couplings of the electromagnetic cur- 
rent to the virtual meson propagating in the intermediate 
meson-baryon state. 

For virtual pion transitions, the LNA terms have the 
very simple form x m ir ~ rriq 1 ' 2 , with values for x as sum- 
marised in Tables IXIVI and [XVI While this contribution 



2. Quark sector magnetic moments 

The u-quark contribution to the proton and E + mag- 
netic moments are illustrated in Fig. 1171 The contri- 
bution Up was described as the most optimal channel 
for directly observing chiral non-analytic curvature in 
quenched lattice QCD simulations [lj| and this curva- 
ture is evident in Fig. El 

The value of x f° r u p is large and negative, predict- 
ing LNA curvature towards positive values as the chiral 
limit is approached. The value of x f° r M s+ vanishes in 
the 7r channel. Similarly, strength in the Sif channel is 
energetically suppressed. Hence the chiral curvature is 
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TABLE XVI: Nucleon magnetic moments and their quark sector contributions. Sector contributions are indicated for single 
quarks having unit charge. Baryon charge states are also summarised. Values for m^ and magnetic moments are in units of 
GeV 2 and /ijv respectively. 

ml u n d„ p n 



0.9972(55) 
0.8947(54) 
0.7931(53) 
0.6910(35) 
0.5925(33) 
0.4854(31) 
0.3795(31) 
0.2839(33) 
0.2153(35) 
0.1384(43) 
0.0939(44) 



0.960(13) 
0.995(15) 
1.032(18) 
1.064(14) 
1.108(17) 
1.163(24) 
1.231(29) 
1.317(41) 
1.395(62) 
1.517(79) 
1.54(11) 



-0.366(8) 
-0.373(9) 
-0.382(10) 
-0.386(10) 
-0.395(12) 
-0.406(16) 
-0.424(20) 
-0.433(28) 
-0.447(47) 
-0.456(60) 
-0.521(87) 



1.401(18) 
1.451(21) 
1.503(25) 
1.547(19) 
1.610(24) 
1.686(33) 
1.783(41) 
1.901(57) 
2.009(87) 
2.17(11) 
2.22(15) 



-0.883(11) 
-0.913(12) 
-0.943(14) 
-0.967(12) 
-1.002(15) 
-1.046(20) 
-1.104(25) 
-1.167(36) 
-1.228(56) 
-1.315(71) 
-1.372(92) 



TABLE XVII: E magnetic moments and their quark sector contributions. Sector contributions are indicated for single quarks 
having unit charge. Baryon charge states are also summarised. Values for m n and magnetic moments are in units of GeV and 
jj,n respectively. 

ml u s+ s E E+ E° E - 



0.9972(55) 


0.976(15) 


-0.416(13) 


1.440(21) 


0.464(7) 


-0.512(10) 


0.8947(54) 


1.010(18) 


-0.414(14) 


1.484(24) 


0.474(8) 


-0.535(12) 


0.7931(53) 


1.044(21) 


-0.412(15) 


1.530(28) 


0.485(9) 


-0.559(14) 


0.6910(35) 


1.072(15) 


-0.408(12) 


1.565(21) 


0.493(7) 


-0.579(11) 


0.5925(33) 


1.113(18) 


-0.407(14) 


1.620(26) 


0.507(8) 


-0.607(13) 


0.4854(31) 


1.163(24) 


-0.406(16) 


1.686(33) 


0.523(10) 


-0.640(16) 


0.3795(31) 


1.221(27) 


-0.409(18) 


1.764(38) 


0.543(12) 


-0.678(18) 


0.2839(33) 


1.286(34) 


-0.412(21) 


1.852(47) 


0.566(14) 


-0.720(22) 


0.2153(35) 


1.344(42) 


-0.421(27) 


1.932(58) 


0.588(18) 


-0.756(28) 


0.1384(43) 


1.418(50) 


-0.421(34) 


2.031(72) 


0.613(23) 


-0.805(32) 


0.0939(44) 


1.446(77) 


-0.429(42) 


2.07(11) 


0.625(30) 


-0.821(53) 



predicted to be negligible for u s + and will contrast the 
upward curvature for u p . This is observed in our lattice 
simulations. Figure 1171 reveals curvature in u p relative 
to a rather linear approach for u s + to the chiral limit. 

The results for u p and u s + are highly correlated and 
therefore the enhancement of the magnetic moment of 
u in the proton over the E + provides significant evi- 
dence of chiral non-analytic behaviour in accord with 
the LNA predictions of chiral perturbation theory. The 
strong correlation of these results is evident in the SU(3) 
flavour-symmetric point at to 2 ~ 0.5 GeV 2 where the 
results are identical. To expose the significance of this 
result, we present Fig. 1181 illustrating the correlated ratio 
of magnetic moment contributions u p /u^+ . There, the 
significance exceeds two standard deviations for quark 
masses between the lightest quark mass considered and 
the SU(3) flavour limit at m 2 ~ 0.5 GeV 2 . 

Figure EI] illustrates the magnetic moment contribu- 
tion of the single u quark in the neutron and the S° , nor- 
malised to unit charge. The magnetic moments match 
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FIG. 18: Ratio of u-quark contributions to the magnetic mo- 
ments of the proton and E + . Chiral curvature in the u-quark 
contribution to the proton's moment gives rise to significant 
enhancement in u„. 
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TABLE XVIII: The A magnetic moment and its quark sector contributions. Sector contributions are indicated for single 
quarks having unit charge. Baryon charge states are also summarised. Values for m„- and magnetic moments are in units of 
GeV 2 and /iat respectively. 



"A 



A'A 



A" 



0.9972(55) 
0.8947(54) 
0.7931(53) 
0.6910(35) 
0.5925(33) 
0.4854(31) 
0.3795(31) 
0.2839(33) 
0.2153(35) 
0.1384(43) 
0.0939(44) 



0.086(10) 
0.091(12) 
0.095(13) 
0.102(11) 
0.109(12) 
0.117(14) 
0.124(16) 
0.135(18) 
0.140(24) 
0.151(27) 
0.154(31) 



1.611(29) 
1.621(20) 
1.631(32) 
1.650(21) 
1.666(24) 
1.688(28) 
1.715(32) 
1.749(37) 
1.780(43) 
1.799(49) 
1.804(53) 



-0.509(9) 
-0.510(10) 
-0.512(10) 
-0.516(8) 
-0.519(8) 
-0.524(10) 
-0.530(11) 
-0.538(13) 
-0.547(16) 
-0.549(18) 
-0.550(20) 



TABLE XIX: S magnetic moments and their quark sector contributions. Sector contributions are indicated for single quarks 
having unit charge. Baryon charge states are also summarised. Values for m 2 and magnetic moments are in units of GeV 2 and 
fj,N respectively. 



0.9972(55) 
0.8947(54) 
0.7931(53) 
0.6910(35) 
0.5925(33) 
0.4854(31) 
0.3795(31) 
0.2839(33) 
0.2153(35) 
0.1384(43) 
0.0939(44) 



1.132(32) 
1.137(33) 
1.141(34) 
1.152(21) 
1.157(22) 
1.163(24) 
1.172(24) 
1.181(25) 
1.189(27) 
1.199(29) 
1.212(29) 



-0.361(13) 
-0.371(14) 
-0.381(15) 
-0.385(13) 
-0.395(14) 
-0.406(16) 
-0.421(18) 
-0.437(20) 
-0.454(23) 
-0.475(27) 
-0.495(31) 



-0.996(23) 
-1.005(24) 
-1.015(24) 
-1.025(18) 
-1.035(19) 
-1.046(20) 
-1.062(21) 
-1.079(23) 
-1.096(25) 
-1.116(29) 
-1.138(32) 



-0.635(22) 
-0.634(23) 
-0.634(24) 
-0.640(14) 
-0.640(15) 
-0.640(16) 
-0.641(17) 
-0.642(17) 
-0.642(18) 
-0.641(20) 
-0.643(20) 
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FIG. 19: The u-quark contribution (single quark of unit 
charge) to the magnetic moments of the neutron, u n , and 
H°, u H o. The magnetic moment for u„ is shifted to the right 
for clarity. 



at the SU (3) -flavour limit where m 2 ~ 0.5 GeV 2 as re- 
quired. The environment sensitivity of the u quark con- 
tribution is subtle and is most evident in the size of the 
statistical error bar. 

The chiral coefficient, %, of the non-analytic term 
~ m^ for u n is large and greater than zero, predicting 
curvature towards negative values as the chiral limit is 
approached. While the coefficient for u s o vanishes in the 
■k channel, a substantial coefficient resides in the energeti- 
cally favoured T,K channel and therefore some curvature 
towards negative values are again predicted as the chi- 
ral limit is approached. Fig. EH is in accord with these 
predictions for the chiral curvature. 

We note that for this case of magnetic moments, the 
anticipated chiral curvature is indeed observed for this 
sector. This contrasts the case of charge distribution 
radii, where chiral curvature was to oppose the Compton- 
broadening of the distribution and was not manifest in 
the simulation results. 

It is interesting to examine the ratio of of singly (u n ) 
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FIG. 20: The ratio of singly (u n ) and doubly (u p ) represented 
quark contributions for single quarks of unit charge to nucleon 
magnetic moments. In the simple SU(6) spin-flavour symmet- 
ric quark model the predicted ratio is constant at —1/2. 



and doubly (u p ) represented quark contributions (for sin- 
le quarks of unit charge) to nucleon magnetic moments 
The SU(6) spin- flavour symmetry of the simple 
quark model provides 
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(5.9) 
(5.10) 



where /^9 M is the constituent quark moment. The quark 
moment pre-factors in Eq. (|5.10(l are respectively, SU(6), 
quark number and charge factors. Discarding quark num- 
ber and charge factors, one arrives at the SU(6) predic- 
tion for u n /u p for single quarks of unit charge of —1/2. 
This prediction is to be compared with Fig. 1201 which 
reveals this ratio to be substantially smaller than the 
SU(6) prediction, even at the SU(3) flavour-symmetric 
limit where m 2 ~ 0.5 GeV 2 . This result is in accord 
with Ref. :4j where this effect was first observed in lat- 
tice QCD. 

The gentle slope of the results in Fig.l^Olat larger quark 
masses suggests that the SU(6) spin-flavour symmetric 
quark model prediction of —1/2 will be realised only at 
much heavier quark masses than those examined here. 

Figure |2"T1 shows the magnetic moment contribution of 
the u-quark sector (or equivalently the <i-quark sector) to 
the A magnetic moment, normalised for a single quark 
of unit charge. In simple quark models, this contribution 
is zero as the u and d quarks are in a spin-0, isospin-0 
state. Our simulation results reveal that the dynamics 
of QCD, even in the quenched approximation, are much 
more complex. The contribution of ua differs from zero 
by more than eight standard deviations at the SU(3) 
flavour-symmetric point, and confirms earlier findings jj| 
of a non-trivial role for the light quark sector in the mag- 
netic moment of A . 
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The chiral coefficient for u\ vanishes in the pion chan- 
nel and has only small strength in the energetically 
favoured NK channel. Hence little curvature is antici- 
pated and this is supported by our findings in Fig. |2*T1 

Turning our attention to the strange quark sectors, 
Figs.l22landl23lpresent results for sa, s-b. and ss magnetic 
moments. In our simulations the strange quark mass is 
held fixed and therefore any variation observed in the 
results is purely environmental in origin. While Ss and 
S£ display only a mild environment sensitivity, sa shows 
a remarkably-robust dependence on its light quark envi- 
ronment. 

Recall that in our examination of the environmental 
flavour-symmetry dependence of the strange quark dis- 
tribution, a strong sensitivity was found. When the en- 
vironmental quarks are in an isospin-0 state in the A, 
the strange quark distribution is broad. On the other 
hand, when the environmental quarks are in an isospin-1 
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state in £ baryons, the distribution radius is significantly 
smaller. It appears that the broad distribution of the 
strange quark in A makes it sensitive to the dynamics of 
its neighbours. 

In the case of strange-quark moments, the LNA con- 
tributions are exclusively from transitions involving the 



kaon. Referring to Table IXTVI one sees that both sa and 
s-s have strong transitions to the energetically favourable 
KN and K £ channels respectively. The coefficients are 
negative such that the virtual transitions will act to pro- 
vide curvature towards positive values, enhancing the 
magnetic moments in these cases. This is not the case 
for S£ where the sign is positive and the transition is to 
the energetically unfavoured K'E, channel. In summary, 
QxPT suggests the magnetic moments for sa and ss will 
display curvature that acts to enhance the magnetic mo- 
ment whereas ss will display little curvature. These pre- 
dictions are exactly as observed in Figs. 1221 and 1251 



3. Baryon magnetic moments 

Figure 1241 depicts the magnetic moments of A , £~ 
and H — . As the magnetic moments of A and S - are 
dominated by the strange quark contribution, these mo- 
ments show only a gentle dependence on the quark mass. 
These contrast E _ where the light d quarks dominate the 
moment. 

The hyperon charge states, S~ and S~, have LNA chi- 
ral coefficients which vanish in quenched QCD. On the 
other hand, the A magnetic moment has some positive 
strength in the energetically favoured NK channel, sug- 
gesting curvature towards negative values as the chiral 
limit is approached. These features are manifest in the 
A and S~ moments of Fig. [2J] where the curvature in 
the A moment towards negative values contrasts the in- 
variance of the S~ moment. The approach of the S _ 
moment to the chiral limit is also fairly linear, in accord 
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FIG. 25: Magnetic moment of E°. 



with a vanishing LNA chiral coefficient. 

Figure [23] presents results for the S° baryon where 
chiral curvature is anticipated to be small. However, a 
comparison of p and E + magnetic moments provides a 
favourable opportunity to observe chiral curvature. The 
proton has a strong negative coupling to the pion chan- 
nel, predicting curvature towards positive values as the 
chiral limit is approached. This contrasts the S + where 
the strong coupling is to the energetically unfavourable 
"E.K channel suggesting a more linear approach to the 
chiral limit. 

Figure I2tj1 depicts the variation of these moments with 
quark mass. These results are highly correlated and 
therefore the enhancement of the magnetic moment of 
the proton over the E + provides significant evidence of 
chiral non-analytic behaviour in accord with the LNA 
predictions of chiral perturbation theory. The strong cor- 
relation of these results is evident in the SU(3) flavour- 
symmetric point at mj ~ 0.5 GeV where the results are 
identical. 
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Figure 1271 reports the magnetic moments of the neu- 
tron and E . The neutron provides a favourable case for 
the observation of chiral curvature associated with the 
pion channel. Similarly the S° has significant strength 
in the energetically favoured T,K channel. In both cases 
the chiral coefficient, \ is positive, predicting curvature 
towards negative values as the chiral limit is approached. 
These predictions are in accord with the observations of 
Fig.E] 

In summary, we have performed an unprecedented ex- 
ploration of the light quark-mass properties of octet- 
baryon magnetic moments in quenched QCD. Figure l2~51 
presents our results in the context of recent state of the 
art results from lattice QCD 0, 0, [IJ ■ The precision af- 
forded by 400, 20 3 x 40 lattices and the efficient access to 
the chiral regime enabled by our use of the FLIC fermion 
action is clear. In every case, the LNA curvature pre- 
dicted by chiral perturbation theory is manifest in our 
results. 
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FIG. 28: Proton magnetic moment in nuclear magnetons is 
compared from a variety of lattice simulations. The solid 
squares indicate our current lattice QCD results with FLIC 
fcrmions. The stars indicate the early lattice results of 
Ref. Q. The crosses (only one point) indicate the results 
of Ref. |jg. The open symbols describe the QCDSF collabora- 
tion results [l2|. Open squares indicate results with (3 — 6.0, 
open triangles indicate those with (3 = 6.2 while the open 
diamonds indicate their results with (3 — 6.4. 



4- Ratio of /j, s - to jUa 

The experimentally measured value of magnetic mo- 
ment of S~ is —0.651 ± 0.003/j.at and that of A is 
—0.613 ± 0.004/ijv, making their ratio greater than 1 at 
1.062 ± 0.012. This has presented a long-standing prob- 
lem to constituent quark models which predict the mag- 
netic moment ratio, S~/A°, to be less than one. 

In the simple SU(6) spin-flavour quark model, the 
magnetic moment of S~ is 



Ms- = gM-S 



1 



(5.11) 



where [i s and fid are the magnetic moments of the con- 
stituent s and d quarks respectively. Since the u-d pair 
in A forms a spin-0 state, the magnetic moment of the A 
has a sole contribution from the s quark 



Taking the ratio yields 





MA 




Us- 


Ms- 

MA 


- = 


4 
3 


1 Md 

3 H s 



(5.12) 



(5.13) 



Now since, the magnetic moment of a charged Dirac par- 
ticle goes inversely as its mass, and since the s and d 
quarks have identical charge, the ratio may be written 



Ms- 
Ma 



1 m 8 

3 m d 



(5.14) 



where m^ and m s are the constituent masses of the d 
and s quarks respectively. Given that m s > m^ it is 
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FIG. 29: Magnetic moment ratio /i s - //ia as calculated in 
quenched QCD. The simple quark model results of 0.836 is 
illustrated by the dashed line while the experimental value of 
the ratio is indicated by the dot-dash line. 

inescapable that this ratio is less than 1 in the simple 
quark model. Indeed, the accepted values of d and s 
constituent quark masses place this ratio at 0.836. 

Figure [55] shows the fi^- //xa ratio as a function of 
quark mass as observed in our quenched lattice calcu- 
lations. Remarkably, the ratio is greater than one at all 
quark masses. 

There are two important aspects of our previous dis- 
cussion that give rise to a result exceeding 1. First, as 
illustrated in Fig. |2J3 we have found that the singly rep- 
resented quarks give a contribution to the magnetic mo- 
ment that is much smaller in magnitude than that of 
the SU(6) quark model prediction. This gives rise to a 
40% reduction in the contribution of the second term of 
Eq. ([5~T4) . 

While this is sufficient to correct the ratio to ~ 1, there 
is a second effect. Namely, the light quark sector makes 
a non-trivial contribution to the A magnetic moment. As 
illustrated in Fig. 1211 this contribution is positive for unit 
charge quarks. Since the net charge of the u-d sector is 
+1/3, the contribution of the s quark in A must have a 
negative value whose magnitude exceeds the observed A 
moment. And this is seen in Fig. [53 There, chiral cur- 
vature in s\ makes the ratio of magnetic moment contri- 
butions sa/se ~ 3/2 as opposed to the SU(6) suggestion 
of 4/3. This resolves the long-standing discrepancy. 



C. Magnetic radii 

Using the values for the magnetic moments obtained 
by scaling the individual quark sector contributions to 
Q 2 = 0, and our values for the form factors at finite Q 2 , 
magnetic radii may be determined in exactly the same 
fashion as the electric radii. 

Analogous to the charge radius, we adopt a dipole form 
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FIG. 30: Magnetic radii of the proton and E + . The latter are 
shifted right for clarity. 



for the Q 2 dependence and define the magnetic radius as 



•it) 



Gm(0) 




Gm(0) 

Gm(Q 2 ) 



(5.15) 



The magnetic radii, (tm)/Gm(0), are tabulated in Ta- 
bic IXXI Figures [201 through [221 display the variation of 
the magnetic radii with m 2 for the octet baryons. 

Figure I5UI depicts the magnetic radii of the proton and 
S + as a function of input quark mass. The somewhat 
subtle differences have a simple explanation in terms of 
the more localised strange quark in S. 

In the proton, the long-range nature of the light-quark 
distributions means that their contributions to the mag- 
netic form factor reduce quickly for increasing momen- 
tum transfers. In the case of S + , which has a broadly dis- 
tributed u quark distribution and a narrowly distributed 
s quark distribution, the reduction in magnitude of the 
form factor is less. Here, the s-quark distribution con- 
tributes positively and remains relatively invariant with 
increased resolution. Thus the S + has a larger form fac- 
tor than the proton at finite Q 2 and hence a smaller mag- 
netic radius. 

Figure l3Tl reports the magnetic radii of the neutron and 
E . Following a similar argument as above, the neutron 
is expected to have a larger magnetic radius than the 5° , 
and this is confirmed in the plot. 

Figure [221 illustrates the magnetic radii of A and E° 
as a function of the input quark mass. In A, most of 
the magnetic moment has its origin in the s quark and 
therefore the magnetic radius will be relatively small. In 
the S° the u-d sector is a major contributor to the form 
factor. As a result, the form factor reduces more at fi- 
nite momentum transfer, which in turn implies that the 
magnetic radius of the S° will be relatively large. 

Figure 13*31 illustrates the magnetic radii of S~, S + and 
E _ . S + is replotted here to facilitate comparison with 
the other two members of the baryon octet. 
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TABLE XX: Magnetic radii {r\ { }/Gm(0) of the octet baryons in fm 2 for various m\ in GeV 2 

p n E+ E° E- A H° 



0.9972(55) 


0.241(7) 


0.240(6) 


0.254(9) 


0.264(9) 


0.236(10) 


0.328(14) 


0.298(13) 


0.337(16) 


0.8947(54) 


0.255(8) 


0.252(7) 


0.265(10) 


0.273(9) 


0.252(10) 


0.328(14) 


0.302(13) 


0.335(17) 


0.7931(53) 


0.269(9) 


0.266(9) 
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0.283(10) 
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0.306(13) 
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FIG. 31: Magnetic radii of the neutron and H°. The latter 
are shifted to the right for clarity. 
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The latter are shifted 



E~ has the largest magnetic radius among the octet 
baryons and this is to be expected based on our consid- 
erations of the origin of the baryon magnetic moment. 
Here, the doubly represented d quark contributes to the 
total baryon form factor with the same sign, whereas the 
strange sector acts to reduce the magnitude of the total 
form factor. Upon increasing the momentum transfer res- 
olution, the ei-sector is reduced dramatically whereas the 
strange sector, acting to reduce the total form factor, is 
relatively preserved, this gives rise to a large drop in the 
total form factor at finite Q 2 and thus a large magnetic 
radius. 



On the other hand, the singly represented d quark 
in the S~ makes only a small contribution to the S~ 
form factor, and therefore the magnetic radius reflects 
the small distribution of the strange quark. 
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VI. SUMMARY 

We have presented an extensive investigation of the 
electromagnetic properties of octet baryons in quenched 
QCD. The development of the O(a)-improved FLIC 
fermion action has been central to enabling this study. 
The FLIC fermion operator is an efficient nearest- 
neighbour fermion operator with excellent scaling prop- 
erties J23 . The vastly improved chiral properties of this 
operator [28| enables the exploration of the electromag- 
netic form factors at quark masses significantly lighter 
than those investigated in the past. The unprecedented 
nature of our quark masses is illustrated in Figs. II Gland 
1281 for the proton charge radii and magnetic moments, 
respectively. 

Central to our discussion of the results is the search for 
evidence of chiral non-analytic behaviour as predicted by 
chiral perturbation theory. We have discovered that all 
baryons having non-vanishing, energetically-favourable 
couplings to virtual meson-baryon transitions tend to be 
broader than those which do not. This qualitative re- 
alisation provides a simple explanation for the patterns 
revealed in our quenched-QCD simulations. 

Of particular interest is the environmental isospin de- 
pendence of the strange quark distributions in A and 
S . When the environmental quarks are in an isospin-0 
state in the A, the strange quark distribution is broad. 
On the other hand, when the environmental quarks are 
in an isospin-1 state in S baryons, the distribution radius 
is significantly smaller. 

Still, evidence of chiral curvature on our large-volume 
lattice is rather subtle in general and absent in the excep- 
tional case of the singly-represented quark in the neutron 
or 5. In this case, the chiral loop effects act to oppose the 
Compton broadening of the distribution. However, it is 
thought that the restriction of momenta to discrete val- 
ues on the finite-volume lattice prevents the build up of 
strength in the loop integral sufficient to counter the nat- 
ural broadening of the distribution as the quark becomes 
light. It will be interesting to explore this quantitatively 
in finite-volume chiral effective field theory. 

In contrast, chiral curvature is evident in the quark- 
sector contributions to baryon magnetic moments. In 
every case, the curvature predicted by chiral perturbation 
theory is manifest in our results. Of particular mention is 
the comparison of the M-quark contribution to the proton 
and S + illustrated in Figs. 1171 and 1181 The environment 
sensitivity of the s quark in A depicted in Fig. [53 is 
particularly robust. 



We find it remarkable that the leading non-analytic 
features of chiral perturbation theory are observed in our 
simulation results. Naively, one might have expected a 
non-trivial role for the higher order terms of the chiral 
expansion which might have acted to hide the leading 
behaviour. However, the smooth and slow variation of 
our simulation results indicate that these higher order 
terms must sum to provide only a small correction to 
the leading behaviour. These observations indicate that 
rcgularisations of chiral effective field theory which re- 
sum the chiral expansion at each order, to ensure that 
higher order terms sum to only small corrections, will be 
effective in performing quantitative extrapolations to the 
physical point. Indeed work in this direction [3, HH Hg 
has been very successful. 

Comparison of our quenched QCD results with experi- 
ment is not as interesting. The chiral physics of quenched 
QCD differs from the correct chiral physics of full QCD 
and our results explore sufficiently light quark masses 
to reveal these discrepancies. The simulation results do 
not agree with experiment, particularly for light quark 
baryons where chiral physics makes significant contribu- 
tions. However, methods have been discovered for quan- 
titatively estimating the corrections to be encountered in 
simulating full QCD and we refer the interested reader 
to Refs. 0j[lii, 38] for further discussion. 

In future simulations it will be interesting to explore 
the utility of boundary conditions which allow access to 
arbitrarily small momentum transfers, providing oppor- 
tunities to map out hadron form factors in detail. Sim- 
ilarly, by calculating near Q 2 = one would have more 
direct access to the magnetic moment. Nevertheless, such 
boundary conditions cannot be seen to substitute for 
larger volume lattices, as the discretisation of the mo- 
menta due to the finite volume of the lattice acts to sup- 
press chiral non-analytic behaviour. Only with increasing 
lattice volumes will the continuous momentum of chiral 
loops be approximated well on the lattice. 
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